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INTRODUCTION 



I . 

Let us begin with a classical example presented in Chap- 
ter 6 of the textbook by Morse and Kimball [Ref. 1], Suppose 
an airplane carries two bombs to attack a railroad track. 

The plane flies along a course perpendicular to the track, 
and drops the bombs . There are three ways to drop two 
bombs: (1) together in salvo, (2) spaced a suitable distance 

apart aiming the midpoint of the pattern at the center of 
the track, and (3) dropping each on a separate run over the 
target . 

We intend to illustrate the method of calculating the 
probability of track destruction in the three cases, the 
purpose being to determine which method of attack is the 
best . 

Let us consider Case (1) , salvo bombing. The two bombs 
leaving the plane simultaneously will hit the ground at some 
distance apart; the impact point being random. For the time 
being the random vector of the impact position from the 
center of impact is termed the ballistic error. The aiming 
is also not free from error. The plane aims at a point on 
the center of the track, but the center of impact will 
deviate from the aim point. This deviation, or aiming error, 
is also random. The impact point of each bomb is, there- 
fore, composed of two errors; first the aiming error which 
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is common to both bombs, and second the ballistic error 
which varies from one bomb to another. 

If the variance of the ballistic error is small while 
that of the aiming error is large, the two bombs will land 
almost at the same point. Therefore it is expected that 
either both will hit or neither will hit the track. To avoid 
the latter situation and to improve the probability of des- 
truction, it is better to spread the landing pattern of the 
bombs. One way is the method of pattern bombing described 
as Case (2) , another is salvo bombing employing bombs 
with a suitable value of variance of the ballistic error. 

Our main concern in this dissertation is maximization of the 
probability of salvo destruction by choosing a suitable 
value of the variance of the ballistic error. In the follow- 
ing we call it the optimal ballistic dispersion. 

The problem described above is not a new one. In an 
interesting talk before a conference at the Ballistic Re- 
search Laboratories held in March 1955 iRef . 2 } , Merritt 
and King stated that an approximate formula for the optimal 
ballistic dispersion was derived by two Englishmen as early 

as in 1936. Following the work in WW II, a number of arti- 

\ 

cles were published on both the salvo and pattern firing 
models, but since the early 1970 's, it seems that these 
models have attracted less attention. Further detail can be 
found in two excellent review papers by Eckler [Ref. 3J and 
Eckler and Burr [Ref. 4] . However, little is known as to the 
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optimization of the ballistic dispersion in salvo models. 

It is the calculation or approximation of the optimal 
ballistic dispersion that is the subject of this dissertation. 

A. SALVO FIRING 

In the following chapters we deal with only two-dimensional 
salvo models. The reason we adopt the two-dimensional model 
is that it is the most frequent case and plays an important 
role in real world applications. One- or three-dimensional 
theory can also be developed to parallel our investigation. 

Suppose that there is a target in a two-dimensional space, 
and a salvo of n weapons is delivered against the target. 
Delivery error relative to the target is assumed to be com- 
posed of two parts, the aiming error and the ballistic error. 

First we aim at the target. With this aiming, the center 
of impact point of n weapons is determined. Let us adopt a 
Cartesian coordinate system Cx,yl such that the origin coin- 
cides with the center of impact. The aiming cannot be per- 
fect. Let the components of the aiming error be -U and -V. 

Then the position of the target with respect to our coordinate 
system is given by (U,V) . It is assumed that the joint 
probability density of U and V exists and is given by fj^(u,v) . 

Now, a salvo of n weapons is fired against the target 
after CU,V) takes a value (u,v) which is unknown to us and 
only is predictable in a probabilistic sense. The components 
of the impact point of the ith weapon are denoted as and 
Y^. It is assumed that (X^,Y^) are independent and identically 
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distributed random variables with joint density f 2 (x,y) . 

In Fig. 1, a typical geometry is shown. The target is at 
(u,v ) , and the impact points of four weapons are scattered 
around the origin. 

As to the target we adopt the so-called point target 
concept; It is assumed that the target is either completely 
destroyed or else undamaged by each weapon. We neglect any 
possible partial damage and its cumulative effect. The 
probability that a weapon landing at (x,y) destroys the tar- 
get at (u,v) is a function of u-x and v-y, denoted as 
D(u-x,v-y) and is called the damage function. 

The probability of destroying the target is calculated 
in two steps. First, suppose that the random variables 
U and V take some values u and v. The conditional proba- 
bility that the ith weapon destroys the target given U = v 
and V = V is 



00 00 

p(u,v) = / / D(u-x, v-y) f 2 (x,y) dxdy . (1.1) 

— oo — oo 

The impact points and are all assumed to be independent, 
so the conditional survival probability of the target given 
U = u and V = V is {1 - p(u,v)}^. Therefore, we obtain the 
probability that the target is destroyed by the salvo of n 
weapons — salvo kill probability — by averaging the conditional 
survival probability with respect to the distribution of u 



9 




X 



Figure 1 
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and V, and then subtracting it from unity; 



oo 



p=l-/ /{I- p(u,v) (u,v) dudv . 



( 1 . 2 ) 



B. UPPER BOUND ON THE SALVO KILL PROBABILITY 

In the example mentioned earlier, Morse and Kimball ob- 
served that the probability of destroying the track by drop- 
ping two bombs in salvo is always smaller than the destruc- 
tion probability of Case (3) , namely the case in which 
the two bombs are dropped independently on separate passes. 
Similar observations are also pointed out in the review paper 
by Eckler and Burr. The observation is valuable, but there 
has been no proof given to this fundamental property dis- 
covered numerically. Indeed, we have the following proposi- 
tion without any specific assumptions on f^(u,v) , f 2 (x,y) 
and D(x,y) . 

Proposition 1.1 

Let the salvo kill probability with n weapons be denoted 



as P 



n 



00 00 



p 



n 



1 - / / {1 - p(u,v) }^f^ (,u, v) dudv . 



(1.3) 



Then, 
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< 1 



(1 - p^) 



n 



(1.4) 



P 



n 



Namely, salvo P^ is always smaller than the kill probability 
of n independent tries. 

Proof 

The proof is immediate from a theorem in the book by 
Hardy and others [Ref. 5] . Define a mean of a non-negative 
function g(x) for a positive r as, 

Mj„(g) = { / g^f , 

where f = f (x) is a weighting function, positive everywhere 
and. 



/ f dx = 1 . 

Theorem 192 from Ref. 5 tells us that if 0 < r < s and 
M (g) is finite, then 



unless g(x) is a constant. 

The weighting function f corresponds to f^(u,v) in our 
proposition, and g corresponds to our l-p(u,v). Therefore, 
for n > 1 , we have 
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Mj^(l-p) < (1-p) 



or 



{M3_(1-p)}'' < (M^d-p)}'' , 



which means 



(1-P,)’^ < 1 - P„ . C 

1 n 

Here is an outline of the succeeding chapters: In Chap- 

ter II, we investigate the simplest salvo model which is 
characterized by circular normal errors and circular Gaussian 
damage function. After that, an approximate formula for the 
optimal ballistic dispersion is introduced in Chapter III; 
we regard the material in Chapter III as the "core" of this 
dissertation. In Chapter IV, the approximate formula is 
applied to the salvo model of Chapter II, and its accuracy 
is studied. In Chapter V, an approximate method of calcu- 
lating the salvo kill probability of the salvo model with 
the general damage function is presented; this case is more 
difficult than the one where the damage function is circular 
normal, so we resort to further approximations. In Chapter 
VI, we investigate the salvo model with the so-called cookie 
cutter damage function, relying heavily on the results of 
Chapters III and V. Finally, salvo model without circular 
symmetry is discussed in Chapter VII. Our concern throughout 
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is the optimization of artificial dispersion. In some cases 
no closed form solution is available, so we offer efficient 
methods for the computation or approximation of salvo kill 
probability with a given dispersion. In other cases we are 
able to offer closed form solutions. 
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II. SIMPLE SALVO MODEL 



In this chapter we deal with the simplest salvo model 
with circular normal errors and a circularly symmetric 
Gaussian damage function. 

Let us adopt a Cartesian coordinate system (x,y) the 
origin of which coincides with the center of impact of the 
n weapons. The target position with respect to this coor- 
dinate system is denoted as U and V in the x and y direction, 
respectively. (U,V) is a random bias common to all the n 
weapons. In this chapter (U,V) is assumed as a circular normal 
variate; namely, we assume that U and V are independent and 

identically distributed normal random variables with mean 0 
2 

and variance a^. The joint density of U and V is, therefore 



The impact point of the ith weapon is denoted as (X^,Y^) . 
It is assumed that (X^,Y^) , i = l,2,...,n are independent and 
identically distributed circular normal random variables 
with joint density function 



(u, v) 





C2.1) 



f 2 (x,y) 




( 2 . 2 ) 
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As to the damage function, we assume that the so-called 
Gaussian damage function with circular symmetry; If the 
target is at (u,v) and a weapon impacts at (x,y) , then the 
destruction probability is given by 

I X 2 - ,2 

(u-x) +(v-y) 

9 2 

D(u-x,v-y) = e . (2.3) 

In the following, let us call this model the Simple Salvo 
Model. 

According to Grubbs [Ref. 6], the Simple Salvo Model was 
first studied as early as in 1953 by H . K. Weiss in BRL Re- 
port No. 879, "Methods for Computing the Effectiveness of 
Area Weapons." The formula for the salvo kill probability 
(2.7) given in the next section is attributed to him. Later, 
Breaux [Ref. 7] found that this formula was not suitable for 
computation when n is large, and gave another method of com- 
putation. However, it seems that there is no published re- 
search investigating this model in full depth; the studies 
usually come to an end when expressions for the kill proba- 
bility are derived. In this chapter, we investigate first 
the property of the kill probability as a function of param- 
eters involved in the model. Methods for computing the kill 
probability and bounds for the kill probability are also 
dealt with in the later sections. 
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A. SALVO KILL PROBABILITY AS A FUNCTION OF n, AND a 

The salvo kill probability is derived in two steps as 
is stated in Chapter I. First the conditional kill proba- 
bility by the ith weapon is calculated given that the random 
bias (U,V) takes a value (u,v) . 

OO CO 

p(u,v) = / / D (u-x, v-y) f 2 (x, y) dxdy . (1.1) 

— 00 —00 

The salvo kill probability is then obtained by averaging 
1 - { l-p(u, v) over the distribution of (U,V) , as 

00 OO 

p = 1 - / / {1-p (u, v) (u, v) dudv . (1.2) 



Using the assumptions (2.2) and (2.3) , we easily calcu- 
late the conditional kill probability (1.1) 



00 00 



p(u,v) 



^ / / 

2ira2 



/ N 2^, ^ 2 2^ 2 

(u-x) +(v-y) X +y 

2 2 
^ 2a 202 



dxdy 



and get 



p (u, v) 




2 ^ 2 
u +v 

2 2 
2(0+02) 



(2.4) 
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Substituting (2.1) and (2.4) into (1.2), the salvo kill 



probability is 



= 1 - 



2vo‘. 



/ / U - 



a 



2 ^ 2 
u +v 

2 2 
2 (a^+ap 



n 



2 ^ 2 
u +v 

o 2 

2a , 



2^ 2 
a +G 2 



dudv 



Conversion from Cartesian coordinates into polar coordinates 
gives 



= 1 - 



1 



I a 

0 



a 



2 2 
2(a^+02) 



n 



2a! 



2 2 
a + 0 ^ 



r dr 



From this we get 



P 




-t/p 




(2.5) 



where 



A = a^/o^, p= (a^ + a^)/a^ . (2.6) 

There are 4 parameters in our model, but the salvo kill 
probability is determined by three, n, o-^/a, and a^/a, or 
n. A, and p as is seen in (2.5) . 

If the integrand in (2.5) is expanded in a binomial 
series and integrated term by term, we obtain the Weiss 



18 



formula for the salvo kill probability. 



P = ^ . ,2.7) 

P is obviously an increasing fxxnction of n, and is de- 
2 2 

creasing in . The latter is easily verified by in- 

2 2 

specting (2.5), where X/p is constant for fixed 02 /^ • The 
2 2 

02/01 -dependence is, however, not obvious. In (2.5), A is 

2 2 

constant for a fixed value of o^/a , and there are two factors 
which involve p in the integrand, with p tending to increase 
P in one and to decrease P in the other. 

2 2 

Table 1 and Table 2 are given here to show the a^/a - and 
2 2 

02 /a -dependence of the salvo kill probability, where n is 

kept constant, n = 2 for Table 1 and n = 16 for Table 2. For 

2 2 2 2 
fixed n and 02 /a , P is decreasing in o^/a as mentioned 

2 2 

above. However, the 02 /oi -dependence is a little different. 

2 2 

It is observed that, generally, P is decreasing in cJ 2 /a for 

2 2 2 2 

small values of o^^/a , but for a larger o^^/a , P increases 
2 2 

when 02 /ol is increased, reaches a maximum value, and then 

decreases. It is also noted that from the viewpoint of real 

2 2 

world applications, the case of large 0 ^/ 0 . with large n is 

especially interesting; e.g., as seen in Table 2, the kill 

2 2 

probability for n = 16 and o-^/o. = 8.0 is improved to 0.5881 

2 2 

when 02 /a is increased to a suitable value. The maximum P 

2 2 

is more than 1.7 times the original value for 02 /a = 0. 
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Table 1. Salvo Kill Probability 

n = 2 



02/a 


0.00 


0.50 


1.00 


2.00 


4.00 


8.00 


0.0 


1.0000 


0.8333 


0.6667 


0.4667 


0.2889 


0.1634 


0.1 


0.9917 


0.8171 


0.6591 


0.4669 


0.2923 


0.1666 


0-.2 


0.9722 


0.7977 


0.6487 


0.4647 


0.2940 


0.1689 


0.3 


0.9467 


0.7767 


0.6365 


0.4609 


0.2946 


0.1706 


0.4 


0 .9184 


0.7550 


0.6232 


0.4560 


0.2944 


0.1717 


0.5 


0.8889 


0.7333 


0.6095 


0.4502 


0.2935 


0.1724 


0.6 


0 . 8594 


0.7120 


0.5956 


0.4439 


0.2920 


0.1728 


0.7 


0.8304 


0.6912 


0.5818 


0.4373 


0.2902 


0.1730 


0.8 


0.8025 


0.6712 


0.5681 


0.4305 


0.2881 


0.1729 


0.9 


0.7756 


0.6518 


0.5547 


0.4236 


0.2858 


0.1726 


1.0 


0.7500 


0.6333 


0.5417 


0.4167 


0.2833 


0.1722 
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Table 2 Salvo Kill Probability 



n = 16 



2 , 2 
Oj/tx 


0.00 


0.50 


1.00 


2.00 


4.00 


8.00 


0 . 00 


1.0000 


0.9935 


0.9412 


0.7835 


0.5511 


0.3370 


0.25 


1 . 0000 


0.9956 


0.9570 


0.8234 


0.6013 


0.3784 


0.50 


1 .0000 


0.9965 


0.9659 


0.8504 


0.6399 


0.4128 


0.75 


1.0000 


0.9968 


0.9710 


0.8692 


0.6700 


0.4417 


1.00 


1 .0000 


0.9967 


0.9739 


0.8824 


0.6938 


0.4662 


1.25 


0.9999 


0.9962 


0.9752 


0.8915 


0.7126 


0.4871 


1.50 


0.9997 


0.9954 


0.9753 


0.8978 


0.7276 


0.5049 


1.75 


0.9993 


0.9941 


0.9745 


0.9017 


0.7394 


0.5201 


2.00 


0.9985 


0.9922 


0.9728 


0.9038 


0.7485 


0.5330 


2.25 


0.9972 


0.9899 


0.9704 


0.9045 


0.7555 


0.5440 


2.50 


0.9954 


0.9869 


0.9672 


0.9038 


0.7606 


0.5534 


2.75 


0.9930 


0.9833 


0.9635 

s 


0.9021 


0.7642 


0.5612 


3.00 


0 . 9900 


0.9792 


0.9591 


0.8996 


0.7665 


0.5678 


3.25 


0.9863 


0.9745 


0.9541 


0.8962 


0.7676 


0.5732 


3.50 


0.9821 


0.9692 


0.9486 


0.8921 


0.7677 


0.5776 


3.75 


0.9772 


0.9635 


0.9427 


0.8875 


0.7670 


0.5810 


4.00 


0.9719 


0.9573 


0.9364 


0.8824 


0.7655 


0.5837 


4.25 


0.9660 


0.9507 


0.9298 


0.8769 


0.7634 


0.5857 


4.50 


0.9597 


0.9438 


0.9228 


0.8710 


0.7607 


0.5870 


4.75 


0.9530 


0.9366 


0.9155 


0.8648 


0.7576 


0.5878 


5.00 


0.9459 


0.9291 


0.9081 


0.8583 


0.7541 


0.5881 


5.25 


0.9386 


0.9214 


0.9004 


0.8516 


0.7502 


0.5880 


5.50 


0.9309 


0.9135 


0.8926 


0.8448 


0.7460 


0.5874 


5.75 


0.9231 


0.9055 


0.8847 


0.8378 


0.7416 


0.5865 


6.00 


0.9151 


0.8973 


0.8767 


0.8307 


0.7369 


0.5853 
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Experience shows that there is one and only one local maximum 

2 2 2 2 

of P as a function of o^/ot for fixed n and a^/a , but we 
have not succeeded in obtaining a proof. 

The salvo kill probability can be increased by increasing 

2 2 2 2 
02/06 from zero for given n and a^/a if the condition below 

is satisfied. 



Proposition 2.1 

The salvo kill probability P for fixed n, and a is 
not maximal at 02 = 0 if and only if 



1 t t 



1+A 2+A 



n+A 



( 2 . 8 ) 



where 



^ 2 , 2 

A = a / 0 ^ . 

Proof 

We will show that 

dP/dp| >0 (2.9) 

p=A 

is equivalent to (2.8) . Using (2.7 ) , (2.9) is written as 

n / (l-e"'*^/^)^“^e“''^/^“^dt < n / (l-e"'*^'^^) ^“^e"^^^“''^t/A dt 
0 0 

( 2 . 10 ) 
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Now define 



I(x) = 



n 



/ 



(1 - e"'*^/^) 



n-l^-t/A-xt 



dt 



Then the inequality (2.10) is rewritten as 

1(1) < -I' (1)/A , (2.11) 

where I' (1) means the value of the derivative dl/dx at x = 1. 
Let s = exp(-t/X) in the integral, then I (x) is 

1 , 

I (x) = nA / s'^^(l-s)^ ds 

0 

_ , n n-1 ^ ^ ^ 1 

Ax+1 Ax+2 Ax+n 



Therefore, 



I' (X) ^ _ y A 

I(x) Ax+j 

The inequality (2.8) is readily obtained from this 
equation with x = 1 and (2.11) . □ 

Proposition 2.2 

A sufficient condition that the salvo kill probability 
P for fixed n, and a is not maximal at C 2 = 0 is 
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n/(e-l) - 1 . 



(2.12) 



2.2 < 

a /a^ ^ 



Proof 

Since l/(x+A) is decreasing in x. 



n 



I 

j=l 



1 

j+A 



> 



n 



I / 

j=i 



j+i 

j 



dx 

x+A 



n+1 

/ 



1 




In 



n+l+A 
1+A • 



If the right-hand-side is equal to or greater than unity, 
i.e., 1+n/ (1+A) = e, then the condition (2.8) in Proposition 
2.1 is satisfied. The relation (2.12) is immediate from the 
above-mentioned inequality. □ 

Proposition 2.3 

Assume that (2.8) holds. A necessary condition that the 
salvo kill probability be maximum at some point <^2 ^ ^ 
fixed n, and a is 

I (-1)^~^(!?) j (~)^~^ Ptj z j- = 0 , (2.13) 

j=l ^ P (p+j) 



where 



, 2,2 , 2 ^ 2 . ,2 

A = a /a-^ , p = (a +a2)/P2_ • 

Proof 

Differentiate (2.7) with respect to p, and equat it to 
0. The equation (2.13) is immediate. □ 
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Example. Case n = 2 

Equation (2.13) is generally too complicated to be solved 
for o^/a in closed form. Here we solve it for the simplest 
case, n = 2. 

The inequality (2.10) with n = 2 is 

1,1 

1+X 2+X ^ ^ ' 

and after a bit of algebra, this reduces to 

1/A = a^/a^ > 2/(/5 - 1) . (2.14) 

2 2 

Now, assiame that o^/a. is large and satisfies (2.14). From 
(2.13) 



P - 2 - 

(p+1)^ ^ (p+2)^ 



0 . 



Therefore, 



2 , 2 

02/a 



= p/A - 1 = 



(p+1) 
p (p + 2) 



- 1 



2 , 2 
01 /a 



= 1/A = 



(p+l)" 

p^(p+2)^ 



(2.15) 



2 2 

This IS a parametric expression of maximizing ^ 2 / 0 . as a 
2 2 

function of 0^^/a the curve of which is shown in Fig. 2. The 
value of the parameter p which corresponds to the point 
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2 / 2 
<r^ / a 



Figure 2. 
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(a^/a^ = {/5+D/2, a^/a^ = 0) is p = Pq 



where 



P 



0 



(/5 - l)/2 . 



When p approaches zero 



2/2 n // 2 

a^/o. ~ 1/4 p 



o^/a^ . ~ 1/4 p , 



and therefore, the asymptotic form of the curve is 





(2.16) 



B. COMPUTING THE SALVO KILL PROBABILITY 

In the preceding section, the Weiss formula for the salvo 
kill probability was given in (2.7). According to Breaxix 
[Ref. 7], this formula is not suitable for calculation for a 
large n because it is an alternating series. Breaux pointed 
out that the salvo kill probability could be expressed in 
terms of an Incomplete Beta Function, and that for large n, 
Tang's method [Ref. 8] would be advisable for calculating an 
Incomplete Beta Function. 

First, we present a recurrence formula for calculating the 
salvo kill probability which is, in principle, the same as 
that of Tang's method. Its generalization will be given in 
a later chapter . 
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Proposition 2.4 

Let us denote the salvo miss probability by n weapons as 
Q(n) = 1-P. 



Q(n) = / e~^ dt . 

0 



(2.17) 



Then Q(n) satisfies the following recurrence equations; 



Q(n) = ^ Q(n-l) , n = 1,2, 



(2.18) 



Q(0) = 1 , 



(2.19) 



where 



U = X/p = 



2 // 2 ^ 2 , 

a / (a + ^2) / P 



/ 2 ^ 2 , , 2 

Ca +Q2)/o-^ 



Proof 

From (2.17) 



Q(n) = / (l-vie~^^^)" e~^ dt 

0 



~t/p,n-l -tj. -t/p.n-^1 -t/p-t,. 

= j (1-ye ' e dt - y j (,1-ye ' e ' dt 

0 0 
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The first term in the right hand side is Q(n-l). In the 
second term, we integrate it by parts, and have 

Q(n) - Q(n-l) - ;^(l-ye “ ^ / (l~ye ^ dt 

0^0 

= Q(n-l) + ^(1-y)^ - ^ Q(n) . 

The equation (2.18) is immediate. D 

It is worthwhile to note that the method presented in 
Proposition 2.4 gives not only the value of Q(n) for a 
specified n^, but also gives all the Q(n) 's up to Uq. It 
is an important advantage of this method over the Weiss 
formula because we are often interested in the salvo kill 
probabilities for several n’s. As is shown in Table 3, the 
computational time for Q(n) by our recurrence algorithm is 
about twice the time required for computation of a single 
Q(n) using the Weiss formula. When the salvo kill probabili- 
ties for more than a single n are needed, therefore, use of 
the recurrence algorithm given in Proposition 2.4 is 
recommended . 

We have another efficient algorithm for computing the 
salvo kill probability, which is good even for large n. 
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Table 3. Comparison of the Processing Time for 
Computing the Salvo Kill Probability 



Mean 


time (microsecond) , 


p = 1 , and y 


= 0.5. 


n 


Weiss formula 


Recurrence 


formula 


10 


92 


175 




20 


189 


383 




30 


300 


587 




Mean 


time (microsecond) , 


p = 1 , and n 


= 10. 




Weiss formula 


Recurrence 


formula 


0.1 


90 


196 




0.5 


92 


175 




0.9 


100 


186 





Computation is on the IBM 3033 installed in the Naval 
Postgraduate School. Each value is the mean of 4 trials. 
One trial is 1000 repetitions of calculation, and the 
elapsed time units are divided by 1000 and multiplied by 
26 (microseconds) . 
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Proposition 2.5 



The salvo miss probability Q(n) by n weapons 



oo 

Q(n) = / (1 - ^ e"'*^ dt 

0 



satisfies the following relationship 



Q(n) = I (J) (1-u)’^"^ q(k) , 

k=0 ^ 

where q(k) is the salvo miss probability by k weapons 
y = 1/ 

OO 

q(k) = / (1 - ^ e"^ dt . 

0 

Proof 

OO 

Q(n) = / {(1-y) + y(l-e P) e ^ dt 

0 



V /n» n . n-k k ,, -t/p.k -t 
I (,) (1-y) y } (1-ye e dt . 

k=0 0 



□ 



(2.17) 



( 2 . 20 ) 



with 
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Corollary 2.6 



The salvo miss probability (2.17) can be calculated by 
the formula 



n 

Q(n) = I A, , (2.21) 

k=0 



where 



= y (n-k+1) / (1-y ) (k+p) , k = 1,2,... 



( 2 . 22 ) 



Aq = (1-y)^ . 



Proof 

As can be seen from (2.18) and (2.19) with y = 1, q(0) 
is unity, and q(k)/q(k-l) is k/(k+p). Utilizing the fact 
that (]^) / = (n-k+1) /k, the corollary can be obtained by 

substitution. □ 

Note that all the terms in the expression (2.24) are 
positive, so there is no problem arising from cancellation 
of terms with alternating signs as was observed in the Weiss 
formula (2.7) . 

In checking the computation, the following P's for special 
cases would be of use. 

(1) The case = 0 corresponds to p -»■ <», and we have 
P = 1 - (1 -y)^ . 
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( 2 ) 



The case = 0 corresponds to y = 1 , p = A , and so 
from the recxirrence formula, we have 



= 1 - 



n(n~l) ♦ 
(n+A) (n-l+A) 



(1+A) 



(2.24) 



2 2 2 

(3) The case a +02 = 0 ^. The formula (2.5) with p = 1 
is easily integrated to give 



P 




(1-A)"^^} 



(2.25) 



where 



y 



2 2 ^ 2 , , 2,2 

a /(a + a^) / A = a • 



C. BOUNDS TO THE SALVO KILL PROBABILITY 

In Chapter I it was pointed out that an upper bound of 
the salvo kill probability is given by the kill probability 
of n repeated independent shots. 



P < P 



where 

C» 00 

P = 1 - [/ / {1 - p (u,v) }f (u,v) dudv]^ . 



Under the assumptions of the present chapter, we have 
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Proposition 2.7 

An upper bound to the salvo kill probability is given by 

P = 1 - (1 - , (2.26) 

where 



2 2 2 9 9 

X = a. /a^ , p = (a • 

As to a lower bound, we have 

Proposition 2.8 (Merritt and King) 

A lower bound to the salvo kill probability is given by 

P = 1 - p(ny)"P ' (2.27) 

where T (p) is the Incomplete Gamma Function, 



r (p) = / e ^ t^ ^ dt , 

^ 0 

and 

2 / / 2 ^ 2, , 2 ^ 2, , 2 
y = a /(a +<^ 2 ^ ' p = (a +o^)/a^ . 

Table 4 and Table 5 show how the true kill probabilities 
are bracketed by these two bounds. The triplets of entries 
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Table 4 Salvo Kill Probability and Its Bounds 

n = 2 



2. 2 
o^/a 


0 . 00 


0.50 


1.00 


2.00 


4.00 


3.00 


0.5 


0.3889 

0.8889 

0.7364 


0.7500 

0.7333 

0.6187 


0 . 6400 
0.6095 
0.5216 


0.4898 

0.4502 

0.3911 


0.3306 

0.2935 

0.2580 


0.1994 

0.1724 

0.1527 


1.0 


0.7500 

0.7500 

0.6321 


0 . 6400 
0.6333 
0.5443 


0.5556 

0.5417 

0.4715 


0.4375 

0.4167 

0.3679 


0.3056 

0.2833 

0.2532 


0.1900 

0.1722 

0.1552 


1.5 


0 . 6400 
0 . 6400 
0.5507 


0.5556 

0.5524 

0.4832 


0.4898 

0.4825 

0.4267 


0.3951 

0.3829 

0.3429 


0.2840 

0.2696 

0.2442 


0.1814 

0.1689 

0.1542 


2.0 


0.5556 

0.5556 

0.4866 


0.4898 

0.4881 

0.4333 


0.4375 

0.4333 

0.3882 


0.3600 

0.3524 

0.3192 


0.2653 

0.2554 

0.2338 


0.1736 

0.1643 

0.1516 


2.5 


0.4898 

0.4898 

0.4353 


0.4375 

0.4365 

0.3922 


0.3951 

0.3925 

0.3554 


0.3306 

0.3255 

0.2977 


0.2489 

0.2418 

0.2233 


0.1664 

0.1593 

0.1482 


3.0 


0.4375 

0.4375 

0.3935 


0.3951 

0.3944 

0.3580 


0.3600 

0.3583 

0.3274 


0.3056 

0.3021 

0.2784 


0.2344 

0.2292 

0.2131 


0.1597 

0.1542 

0.1444 


3.5 


0.3951 

0.3951 

0.3588 


0.3600 

0.3596 

0.3291 


0.3306 

0.3294 

0.3032 


0.2840 

0.2815 

0.2611 


0.2215 

0.2175 

0.2034 


0.1536 

0.1492 

0.1404 


4.0 


0.3600 

0.3600 

0.3297 


0.3306 

0.3303 

0.3045 


0.3056 

0.3048 

0.2823 


0.2653 

0.2635 

0.2457 


0.2099 

0.2068 

0.1944 


0.1479 

0.1443 

0.1365 



Triplets are, from the top, P, the upper bound, 
P, the true value, and P, the lower bound. 
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Table 5 Salvo Kill Probability and Its Bounds 



n = 16 



2/ 2 
Oj/a 


0.00 


0.50 


1.00 


2.00 


4.00 


8.00 


1.0 


1 . 0000 
1.0000 
0.9997 


0.9997 

0.9967 

0.9944 


0.9985 

0.9739 

0.9688 


0 . 9900 
0.8824 
0.3750 


0.9459 

0.6938 

0.6867 


0.8147 

0.4662 

0.4611 


2.0 


0.9985 

0.9885 

0.9952 


0.9954 

0.9922 

0.9862 


0 . 9900 
0.9728 
0.9644 


0.9719 

0.9038 

0.8935 


0.9151 

0.7485 

0.7388 


0.7824 

0.5330 

0.5258 


3.0 


0.9900 

0.9900 

0.9817 


0.9821 

0.9792 

0.9685 


0.9719 

0.9591 

0.9469 


0.9459 

0.8996 

0.8864 


0.8819 

0.7665 

0.7546 


0.7515 

0.5678 

0.5590 


4.0 


0.9719 

0.9719 

0.9592 


0.9597 

0.9573 

0.9432 


0.9459 

0.9364 

0.9215 


0.9151 

0.8824 

0.8674 


0.8481 

0.7655 

0.7522 


0.7222 

0.5837 

0.5739 


5.0 


0.9459 

0.9459 

0.9305 


0.9309 

0.9291 

0.9130 


0.9151 

0.9081 

0.8917 


0.8819 

0.8583 

0.8424 


0.8147 

0.7541 

0.7402 


0.6945 

0.5881 

0.5778 


6.0 


0.9151 

0.9151 

0.8983 


0.8987 

0.8973 

0.8802 


0.8819 

0.8767 

0.8597 


0.8481 

0.8307 

0.8145 


0.7824 

0.7369 

0.7229 


0.6684 

0.5853 

0.5749 


7.0 


0.8819 

0.8819 

0.8647 


0.8650 

0.8639 

0.8468 


0.8481 

0.8441 

0.8273 


0.8147 

0.8016 

0.7857 


0.7515 

0.7167 

0.7030 


0.6439 

0.5780 

0.5677 


8.0 


0.8481 

0.8481 

0.8310 


0.8313 

0.8305 

0.8136 


0.8147 

0.8117 

0.7953 


0.7824 

0.7724 

0.7570 


0.7222 

0.6952 

0.6819 


0.6209 
0 . 5680 
0.5579 



Triplets are, from the top, P, the upper bound, 
P, the true value, and P, the lower bound. 
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in the table are, from the above, the upper bound P, true 

value P, and P, the lower bound. It is observed that 
— 2 2 

P - P is not small for smaller values of ^2^^ • gets smaller 

2 2 2 2 — 
when is increased. If 0^/0. is small, P and P are 

close. On the other hand, P is fairly close to P over wide 
2 2 2 2 

ranges of o^/ot and ^2/01 . Figures 3 and 4 illustrate this. 

It seems reasonable that Merritt and King used P as an 
approximate formula for P in their study on the optimal 
ballistic dispersion, because computation was really a 
problem at that time. Nowadays, however, direct computation 
of P is much more natural. 
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<r /a 
2 

Figure 3 '. 
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^2/ a 2 



Figure 4. 
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III. OPTIMIZATION OF SALVO FIRING 



In Chapter II we have learned the following: The salvo 

kill probability is a decreasing function of if other 
parameters are kept constant. Therefore, to obtain higher 
salvo kill probability, it is necessary to make as small 
as possible. The salvo kill probability as a function of 
02 is, however, complicated. If is sufficiently small, 
the salvo kill probability is monotone decreasing in 02/ 
but if 0j^ is not so small and satisfies the condition (. 2 . 8 ) , 
then the salvo kill probability is not maximal at 02 = 0/ 
but takes a maximum value at some 02 > 0 . In other words, the 
smallest ballistic dispersion is not always most desirable. 
Deliberate increase of the ballistic dispersion is encoun- 
tered in shotguns and in related military weapons. In this 
chapter we present an approximate formula for the optimal 
ballistic dispersion. We will first briefly sketch Walsh's 
theory [Ref. 9 ] since our approximate formula is based on 
it. In Section B our approximate formula is presented. 

Its accuracy will be studied in Chapter IV. 

A . WALSH THEORY 

The Walsh model is based on assumptions similar to those 
of Chapter I. A salvo of n weapons is fired at a target. 

The impact point error is composed of two parts, a random 
bias which is common to all n weapons, and a round- to-round 
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■mm 



M 




m: 



error. A Cartesian coordinate system is chosen so that its 
origin coincides with the center of impact of n weapons. 

The coordinates of the random position of the target are 
denoted as U and V with respect to this system. The joint 
probability density of U and V is given by f^(u,v) . The 
impact points of weapons (X^,Y^) are independent and iden- 
tically distributed random variables with joint density 
f 2 (x,y). It is assumed that impact of a weapon at (x,y) , 
given the target is at (u,v) , destroys the target with 
probability D(u-x,v-y) independently of other weapons. 

The lethal area of a weapon is denoted as A: 



00 oo 

j j D(x,y)dxdy = A . 

— OO — oo 



The salvo kill probability is given by (1.2), 



00 00 

P = l- / / {1- p(u,v) }^f^(u,v)dudv , (.1.2) 

— OO —00 



where p(u,v) is the conditional destruction probability by 
a single weapon given that the target is at (u,v) . 



OO 00 

p(u,v) = / / D (u-x, v-y) f 2 (x,y) dxdy . (1.1) 

— OO — oo 



Note that p(u,v) can be integrated to give A. 
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A . 



(3.1) 



OO OO 00 00 

/ / p(u,v)dudv = A / / f 2 (x,y)dxdy = 

— OO —00 — OO —00 



The problem is to derive the optimal ballistic disper- 
2 

Sion, 02* Instead of solving the original problem, Walsh 

solved a revised problem that was different from the original 

2 

one in that the variable is not 02 itself but rather the 
function p(u,v) . 

Proposition 3.1 (Walsh) 



00 00 



max P 
p( • ♦ ) 



= 1 - / / {1 - p(u,v) }^f (u, v) dudv (3.2) 



— OO —00 



subject to 



00 00 



/ / p(u,v)dudv = A , 



—00 —00 



(3.3) 



p(u,v) = 0 



has a solution 



1 - {-E — 7 ^ ^ (u,v) € E 

f (u, v) 



p*(u,v) 



C3.4) 



(U,v) € E 



where 



42 



E 



(3.5) 



= { (u,v) 1 f (u, v) = C} . 

C is positive and is given by 

/ / p*{u,v) du dv = A . (3.6) 

E 

The corresponding salvo kill probability is 

P* = CA + / / {f, (u,v) - C} dudv . (3.7) 

E 

In case the random target bias is given by a normal 
distribution, we also have 

Corollary 3.2 

OO 00 

max P = l- J / {1- p (u, v) }^f (u, v) dudv (3.8) 

p ( • • ) —00 — OO 

subject to 

00 OO 

/ I p{u,v) dudv = A , 

—00 — OO 

(3.9) 

p(u,v) ^ 0 , 



where 



u 



2a 



(u, v) 



u 



2a' 



V 



2^ '^u'^v 



(3.10) 
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has the solution 



2 2 



1 - exp{-(j) 4- Yin-l) (u,v) e E 



p*(u,v) 



a o 
u V 



(3.11) 



where 



(u,v) e E 



E = {(u,v) |u^/a^ + v^/a^ 1 2(n-l)()>} 



(}) is positive and is determined by 



e - 1 + 4> - A/2 tt (n-1) a o = 0 . 

^ ' u V 



(3.12) 



The maximim value of the salvo kill probability is given by 



= 1 - '^{1 + (n-1) (1 - e"*^) } . 



(3.13) 



Proof 

The equation (3.11) is readily obtained from (3.4) in 
Proposition 3.1 with f^(u,v) given by (3.10) and C, 



1 

2t\ o a 
u V 



(3.14) 



Substituting (3.11) in (3.6) and changing variables into 
(r,9) by 
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u = /2(n-l) a^r cos 6 , v = /2(n-l) r sin 0 (3.15) 



we have 



A = / / p*(u,v) du dv 

E 

, , 2 

= 2TT(n-l)a a / (1 - e"*^ ^ ) 2r dr 

u V 

= 2ir (n-1) (<}) - 1 + e"**^) . 



Substituting (3.10) and (3.14) in (3.7) with (3.15), 






V * -- 



2ira 



A -(n-l)(|) , 2 tt , 1 ^ , ,, 2 , ,, , 



U V 



U V 



= 1 - + (n-l)(}) - A/27ra a } . 



u V 



This formula together with (3.12) gives (3.13). □ 

For simplicity, let us call the function p*(u,v) the 
(Walsh) optimal p function, and the associated salvo kill 
probability P* the Walsh optimal P. It should be noted that 
the constraint in the revised problem is weaker than that of 
the original problem, and so the Walsh optimal P is an upper 
bound to the maximum salvo kill probability of the original 
problem. It should also be noted that the Walsh optimal p 
function is not feasible in the original problem (being non 
zero on only a finite set) if dispersion is present. 



) 2rdr 
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B. APPROXIMATE FORMULA FOR THE OPTIMAL BALLISTIC DISPERSION 
Walsh developed the theory stated above, but he did not 
extend his theory beyond the results (3.4)-(3.7). The impor 
tant part left by Walsh and subsequent workers is to bridge 
the gap between the revised problem and the original prob- 
lem. Walsh got the optimal p function as given in (3.4), 
whereas we need the optimal density function f 2 (x,y). The 
material in this section bridges that gap. 

Let us define three double Fourier transforms as follows 

oo 00 iCuj.u+w-v) 

p*(tUj^,aj 2 ) = / / ® p*(u,v)dudv, 

— OO — oo 

~ oo oo i(o)^x+uj2y) 

D(uj^,a) 2 ) = / / ® D(x,y)dxdy, (3.16) 

— OO — oo 

~ oo oo i (oj^x+a)2y) 

= / / e ^ f 2 (x,y)dxdy . 

— OO — oo 

Then from (1.1) we have 

p*(w^,uj2) = D(uj^,ai2) £2 '^^2^ * (3.17) 

In the Walsh theory, D(x,y) is a given function, and 
the Walsh optimal p function is given by Proposition 3.1. 
Therefore, both D(w^,W 2 ) and p*(a)^,(iJ 2 ^ (3.17) are now 

known functions. Thus the Fourier transform of the "optimal 
f^(x,y) can be calculated by 



46 



±2 



p* (o)^,co2) 



(3.18) 



D(o)^,aj2) 

and we will get f 2 (x,y) as its inverse Fourier transform. 

This might have been the end of Walsh theory, but there are 
difficulties. First, we do not have any guarantee that 
f 2 (x,y) is a density function, because Walsh solved the 
restricted optimization problem subject to p(u,v) ^ 0 but not 
f 2 (x,y) ^ 0. Further, even if this problem were solved, there 
would still remain the problem of feasibility of the f 2 func- 
tion. Walsh suggests in his paper that it would be possible 
to obtain a real-world dispersion function f 2 (x,y) as the 
theory requires, at least approximately, by careful design 
of the ammunition, but this may be asking too much of the 

manufacturing process. We assume below that the function 

2 

f 2 (x,y) is normal, with only the variance O 2 subject to 
manufacturing control . 

The problem, then, is to determine or approximate the 

2 

optimal ballistic dispersion 02 , given the damage function 
D(x,y) and the function p*(x,y) from (say) Proposition 3.1. 

In the following, both D(x,y) and p*(x,y) are assumed to be 
even functions of x and y respectively, and to have "moments" 
of all degrees. The lethal area is denoted as A. 



00 00 00 00 

/ / xD(x,y)dxdy = / / yD(x,y)dxdy =***=0, (3.19) 

—00 — »oo —00 —00 
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00 



00 



00 



/ / xp*Cx,y)dxdy = j f yp*(x,y)dxdy = ••• = 0 , 



— OO —00 



— oo — oo 



00 oo 



/ / D(x,y) dxdy = A , 



— OO — oo 



00 oo 



/ / x^^D(x,y) dxdy < <», k = 1,2, 



—00 — oo 



00 00 



/ / x^^p*(x,y) dxdy < <», k = l,2, ... 



(3.20) 



—00 — oo 



etc . 



For the sake of generality, we actually consider the non- 
symmetric case where the optimal dispersion may be differ- 
ent in the x and y directions. Let the two dispersions be 

0^ and a^. Then 
X y 



f2 (X/y) 






2 2 
-JL- 
- 2 „ 2 
2a 2a 
. X y 



(3.21) 



with the double Fourier transform 



f 2 (Wj^, (JO 2 ) 



= e 



^22 „2 2 

a u), a 10 - 
X 1 y 2 

■ 2 ■ 2 



( 3 . 22 ) 
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Expanding the exponential function in the double Fourier 



transform of D(x,y), and integrating term by term, we have 



i(to,x+co-y) 

= / / e D(x,y)dxdy 



—CO — oo 



oo . J oo oo 

= 1^1 I (C0j^x+O)2y) (x,y) dx dy 

j = 0 ^ 



= A[l--^|-/ / x^D(x,y)dxdy 



2 

;i 

2 A 
2 

2 A 



— OO —CO 



T i f / y^D ( X , y) dx dy + . . . ] 



— OO — oo 



where the third line is obtained because D is even. There- 

2 2 

fore, for sufficiently small values of and 0 ^ 2 / 






(3.23) 



where 



Dq ((0j^,u)2) 



A exp t — ^ i / / x^D(x,y) dx dy 



— OO —00 



2 
|2 
2 A 



^ i- / / y^D(x,y) dxdyj 



— OO —00 



(3.24) 



Similarly, 
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(3.25) 



p*(oa^,oj2) = Pq(oj^,W2^ 



2 2 

for small oj^ and where 



‘"l 1 r“ r 2 



Pq(oj^,oj 2 ) = A exp t~ ~ ^ / / ^ p*(x,y) dx dy 



—00 — oo 



<"2 1 r 2 



- — ^ / / y P*(x,y) dxdy] 



— 00 — oo 



(3.26) 



2 2 

Now we determine two parameters a and a in (3.21) by 

X y 

the relation 



f2(w^,W2) 



Pq 

Dq 



(3.27) 



Then from (3.22), (3.24), and (3.26) we get an approximate 

formula for the optimal ballistic dispersion; 



00 00 

= [^ / / x^{p* (x,y) -D(x,y) } dx dy, 0]'*’ 

— OO — <x> 



(3.28) 



OO 00 

Oy = [^ / / y^{p* (x,y) -D(x,y) } dx dy, 0]'*' 



—00 — OO 



where 



[a,0] 



. ( 



if a ^ 0 



otherwise 



(3.29) 
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Note 1 



In this proposition, the approximate formula for the 
optimal dispersion is derived from the condition that 

coincides with f2(o)^,a)2) fot infinitesimally small 

2 2 

and 0)2* It would be better if we could derive the approxi 

mate formula from the condition that f 2 (03^,032) is nearly 

2 2 

equal to over a wider domain of and 0)2 which 

include ( 0 , 0 ) . In the next chapter we will give an idea of 
improving the approximate formula taking this into account. 



Note 2 



Generally, approximation ( 3 . 23 ) with ( 3 . 24 ) is accurate 
only up to the second power of o)^ and 002* If the damage 
function is of Gaussian type. 



D(x,y) 



2 2 

^ ^ 

o 2 » 2 

2cc 2a 



r 



then 



DQ(a)j^,o) 2 ) = 0(03^^,032) 

is an exact relationship. 

As a corollary, we have the following approximate formula 
for the model with normal errors and Gaussian damage function 

Consider a salvo firing of n weapons. The density func- 
tions of the random bias and the round-to-round ballistic 
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error are, respectively. 



u 



f^(u,v) 



2a 2a' 

u ■' 



2 IT a a 

u V 



and 



X 



f2(x,y) 



2 IT a a 
X y 



2a 2a 
X y 



the damage function is given by 



(u-x) _ (v-y) 



2a' 



D(u-x,v-y) = e 



u 



2o‘ 



(3.32) 



(3.33) 



(3.34) 



Then approximate formulae for the optimal ballistic dispersion 

for fixed n, a , a , a , and a are 
u V X y 

(a^/ax)^ = [i(n-l)^(t.^(ayax)^(4)ya^)-(n-l) (aya^)^-l,0]'^ 

(3.35), 

(ty/a^)^ = [j(n-D ^4)^ (a^/Uy) ^ (a^/a^) - (n-1) (a^/tty) ^-1, OJ'*’ 

where <{> is determined by the equation 

- I + - (a /a ) (a /a )/(n-l) = 0 . (3.36) 

^ x u y V 
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An upper bound to the maximum salvo kill probability is given 
by 



P = 1 - + (n-1) (1 - e"'^) } . (3.37) 

The solution to the Walsh problem with the above men- 
tioned assumptions is already given in Corollary 3.2. By 
use of (3.11) , we have 

, 00 

^ / / up* (u, v) du dv 

— oo — oo 



3 

a a 

u V 



/2(n-l) 4> T . -<t> + Yll — iY 

2Trj ^ril-e } r 



dr 



3 

TTO a 

U V I 



2,2 



-{ (n-1) <t> - (n-l)A/ira a } 

A Ll V 



and 



i / / x^D(x,y) dxdy = . 

—00 —00 



Substituting these relations into (3.28), we get the approxi- 

2 

mate formula for (3.35). Similar arguments follow for 



53 



IV. SIMPLE SALVO MODEL REVISITED 



In this chapter we apply the general approximate formula 

for the optimal ballistic dispersion obtained in the last 

chapter to the Simple Salvo Model dealt with in Chapter II 

and investigate its accuracy. 

The assumptions in this chapter are the same as those in 

Chapter II. We consider a salvo of n weapons against a 

target. The impact point error is composed of two parts, 

the random bias and the round-to-round error. The position 

of the target (U,V) has a circular normal distribution centered 

2 

at (0,0) with variance . The impact points of the ith 
weapon (X,Y) are assumed independent and identically dis- 
tributed circular normal random variables centered at (0,0). 

2 

with variance 02 * Finally, we assume the circular Gaussian 
damage function with parameter : 




1 




(4.1) 




/ 



f2 (x,y) 



1 




(4.2) 




/ 



(u-x) (v-y) ^ 



D (u-x, v-y) 



e 



(4.3) 
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In (3. 35) , let 



u 



= a = 

V 1 



a = a = 

X y 2 



a 



X 



= a = a 

y 



(4.4) 



then we have the following: 

An approximate formula for the optimal ballistic disper- 
sion of the Simple Salvo Model is given by 



2/2 r ^2 2 , - , + 

a 2 /a = l<p z - z - 1, 0] 



(4.5) 



where 

z = (n-l)a^/a^ = 1/ (e"'^’ - 1 + cj)) . (4.6) 

An upper bound to the maximum salvo kill probability is given 
by 



P = 1 - + (n-1) (1 - e“^) } . (4.7) 

2 2 2 2 
The true optimal O 2 / 0 . is a function of n and o ^/a , but 

2 2 

according to (4.5) and (4.6) the approximate ^ 2 ^'^ deter- 

2 2 . 

mined by z = (n-Da^^/a , (p being a parameter. To obtain an 

2 2 

explicit formula for 02 /a in z instead of the parametric 
form, we have 
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Leirana 4 . 1 



The series solution to the equation in <t> , 

e"'^ - 1 + 4) - u^/2 = 0 (4.8) 

is given by 

= u^ + u^/3 + u'^/12 + u^/60 + ... (4.9) 

Proof 

Let the left-hand-side of (4.8) be expanded in <)) . 

Then 

u^ = 4)^(1 - 4)/3 + 4)^/12 - 4)^60 + ...) . 

If the terms in the right-hand-side of (4.9) are expanded 
in 4> f we have 



2 .2 ^3.. 

u = 4) - 4> /3 


+ 4> /i2 — (p /60 + . . . 


3 3 

u-^/3 = 4)-^/3 


- 4)^/6 + 4>^/18 - . . . 


4 

u /12 = 


(p^/12 - 4>^/18 + . . . 


u^/60 = 


4>^/60 - . . . 



from which (4.9) is obvious. □ 
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Using the Lemma and the relationship (4.6), we can express 

2 2 r—t— 

(4.5) as a power series in /2/z. If we keep only 
the first three terms in the expansion, we get an alterna- 
tive approximate formula for the optimal ballistic disper- 

2 

sion which we will denote as a . 

a 

= I /i72 - I + A /m , (4.10) 



with 



2 = (n-l)a^/a^ . 



In Table 6 two values of the approximate optimal ballis- 
tic dispersion are compared. One is calculated with (4.5) 
and (4.6), the other with (4.10). Agreement seems quite 
satisfactory, differences appearing at most in the fourth 
decimal digit. The formula (4.10) needs more terms to be a 



Table 6. Comparison of Two Formulae 



/ in 2,2 2,2 „ 2.2 

2 = (n-DUj^/a ^ 2 ^ 





800.0 


12.5017 


12.5017 




400.0 


8.5971 


8.5971 




200.0 


5.8367 


5.8367 




100.0 


3.8855 


3.8854 




50.0 


2.5069 


2.5067 




25.0 


1.5335 


1.5331 




12.5 


0.8475 


0.8467 


2 / 2 
02/a 


is calculated by 


(4.5), (4.6), 


and 


,,2/ 2 
Ca/a 


is calculated by 


(4.10) . 
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good approximation for smaller values of z, but later we 
will see that the case with z less than 25 is not important 
from the viewpoint of application. 



Note 



For n = 2, the formula (4.10) gives 




'T 5/2 

T " 6 lO" “'^^l * 



(4.11) 



In Chapter II, the exact formula of the optimal ballistic 
dispersion was derived for n = 2. The series expansion of 



It is interesting to note that every term in (4.11) is smaller 
than the corresponding term in (4.12) : The approximate 

value given by (4.11) for n = 2 is always smaller than the 
true value. 

For comparison's sake, we will also consider the classi- 
cal approximate formula by Merritt and King [Ref. 2] who 

2 

approximate the optimal ballistic dispersion by cr^ where 



2 , 2 . 

IS 





(4.12) 





(4.13) 



z 



c 




(4.14) 
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A. ACCURACY INVESTIGATION 



In this section we will investigate the accuracy of both 

our approximate formula and the classical one. 

Table 7A for z = 400, Table 7B for z = 100, and Table 

7C for z = 25 illustrate the calculation carried out for the 

accuracy study. For a fixed z, the approximate value for 

2 

the optimal ballistic dispersion divided by a is a constant 

which is given next to the value of z in the table. On the 

other hand, the true optimal value increases with n: The 

first column is n and in the second column is given the corres- 

2 2 

ponding true optimal 02 /d , the calculation of which is based 

on the necessary condition of optimality (2.13), and requires 

fairly long calculation time. The relative error of our 

2 2 2 2 

approximate formula Ac = (a x 100% is given in the 

3rd column. Its value is found to always be negative. The 

salvo kill probability associated with = 0 is denoted as 

Pq and is in the fourth column. The fifth column is the upper 

bound P of the maximum salvo kill probability given by (4.6). 

and (4.7) . In the sixth column, we have the true maximum 

salvo kill probability designated simply as P which is calcu- 

2 

lated by (2.21), (2.22) with the optimum 02 given in the 

2nd column. The entries in the 7th column are the salvo kill 

probabilities associated with the approximately optimal dis- 
2 

persion a . The notation P is used for it. The last column 

3 . 3 

is the relative loss in the kill probability which would be 

2 

caused if we were to use the approximate value a instead of 



59 



Table 7 Accuracy of Approximation (4.10) 



A z = 400, 



n 






Po 


5 


9.34 


-7.9 


0.0225 


0 


9.46 


-9.1 


0.0634 


5 


9.57 


-10.2 


0.1088 


0 


9.67 


-11.1 


0.1556 


5 


9.76 


-11.9 


0.2024 


0 


9.84 


-12.6 


0.2484 


5 


9.92 


-13.3 


0.2931 


0 


9.98 


-13.9 


0.3362 



B z = 100 



n 






Po 


4 


4.31 


-9.9 


0.0600 


8 


4.41 


-11.8 


0.1703 


2 


4.48 


-13.4 


0.2828 


6 


4.55 


-14.6 


0.3877 


:0 


4.61 


-15.7 


0.4817 


14 


4.66 


— 16 . 6 


0.5641 


!8 


4.70 


-17.4 


0.6352 


12 


4.74 


-18.1 


0.6959 



C z = 25, 



n 




^ % 


Po 


2 


1.77 


-13.6 


0.0573 


4 


1.82 


-15.9 


0.2137 


6 


1.86 


-17.8 


0.3706 


8 


1.90 


-19.3 


0.5076 


L0 


1.93 


-20.5 


0.6206 


L2 


1.95 


-21.5 


0.7107 


14 


1.98 


-22.4 


0.7812 


L6 


1.99 


-23.1 


0.8356 



= 8 . 60 



P 


P 


Pc. 




0414 


0.0410 


0.0410 


0.0 


1484 


0.1456 


0.1454 


0.1 


2777 


0.2709 


0.2704 


0.2 


4063 


0.3952 


0.3943 


0.2 


5229 


0.5081 


0.5067 


0.3 


6230 


0.6056 


0.6038 


0.3 


7061 


0.6873 


0.6851 


0.3 


7733 


0.7542 


0.7518 


0.3 



= 


3.89 






P 


P 


P« 


^?*/„ 


0905 


0.0891 


0.0890 


0.1 


2947 


0.2872 


0.2866 


0.2 


4952 


0.4809 


0.4794 


0.3 


6558 


0.6371 


0.6350 


0.3 


7727 


0.7529 


- 0.7503 


0.3 


8534 


0.8348 


0.8321 


0.3 


9071 


0.8910 


0.8885 


0.3 


9419 


0.9288 


0.9266 


0.2 



= 


1.53 






P 


P 


P^ 




0660 


0.0652 


0.0651 


0.1 


2733 


0.2663 


0.2657 


0.2 


4822 


0.4674 


0.4658 


0.3 


6497 


0.6299 


0.6274 


0.4 


7710 


0.7497 


0.7468 


0.4 


8539 


0.8339 


0.8309 


0.4 


9085 


0.8912 


0.8884 


0.3 


9435 


0.9295 


0.9271 


0.3 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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2 

the true optimal o , AP = (P-P )/pxioo%. It is noted that 
in the case with a large z value the salvo kill probability 
is much improved by employing the optimal ballistic disper- 
sion. For instance, when z = 400, the ratio P/Pq is around 
2.5. But the improvement is not so much when z gets small. 
In the case with z = 25, the ratio is at most 1.25. 

As to our approximation, we observe the following. 

2 

(1) The approximate a is always smaller than the true 

3 

value, at least for the values of z investigated. 

(2) The discrepancy between the approximate and the 
true values is larger for smaller z. 

2 

(3) For a fixed z, the optimal ballistic dispersion a 2 
increases as n is increased. Thus the discrepancy 
grows when n increases for a given z . 

(4) However, we may say that the discrepancy is not so 
serious because the kill probability associated with 
the approximate optimal ballistic dispersion is not 
so different from the maximxxm value. The relative 
loss given in the last column is at most 0.4% which 
occurs when z is 25, a less important case. 

Figure 5 is presented here to illustrate the properties 

2 

1 and 2. The ordinate is the relative error Aa , but it 

3 

2 

needs several words. The relative error Aa varies with n 

3 

for a fixed z as seen in Table 7. To get a representative 

2 

Aa^ value for a z, an n is picked for the given z such that 
the maximum salvo kill probability is close to 0.6, an 
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Figure §. 
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2 

arbitrarily chosen figure. Then the corresponding Aa is 

Si 

chosen for that z value, for instance -12.63 for z = 400, 
and similar procedure follows for other z's. The bottom 
curve in Figure 5 shows this relationship. Later we will 
discuss other curves in the same figure. 

From the figure we see that the relative error is always 
negative and its absolute value is decreasing in z. The 
z-dependence is, however, relatively small, and we learn 
its bias to the negative side is a characteristic of this 
approximate formula. 

Similar investigation was carried out on the classical 

approximate formula (4.13). Table 8A for z^ = 400, Table 8B 

for z = 100, and Table 8C for z =25 correspond to Table 7, 
c c 

the only difference being that z^ is tabulated, rather than 

z. From the tables we observe the following. 

2 2 

(1) The approximate o^/a is larger than the true value 
except for the cases of very large z^'s. 

(2) The discrepancy between the approximate and the true 
value is larger for smaller z^. This tendency is 
similar to that of our approximation. 

2 

(3) For a fixed z^, the optimal ballistic dispersion 02 
increases as n is increased. Thus, the discrepancy 
gets smaller when n is increased for a given z^ less 
than 800. 

(4) However, the discrepancy is not serious because the 
kill probability associated with the approximate 
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Table 8 Accuracy of Approximation (4.13) 





A 


z ^= 400, 




10.00 










Po 


P 


P 


Pc 


4 P % 


9.34 


7.1 


0.0225 


0.0414 


0.0410 


0.0410 


0.0 


9.46 


5.7 


0.0634 


0.1484 


0.1456 


0.1455 


0.0 


9.57 


4.5 


0.1088 


0.2777 


0.2709 


0.2709 


0.0 


9.67 


3.4 


0.1556 


0.4063 


0.3952 


0.3951 


0.0 


9.76 


2.5 


0.2024 


0.5229 


0.5081 


0.5080 


0.0 


9.84 


1.6 


0.2484 


0.6230 


0.6056 


0.6056 


0.0 


9.92 


0.9 


0.2931 


0.7061 


0.6873 


0.6873 


0.0 


9.98 


0.2 


0.3362 


0.7733 


0.7542 


0.7542 


0.0 





B 


z ^= 100 




5.00 








% 


Po 


p 


P 


Pc 


^ P/o 


4.31 


15.9 


0 . 0600 


0.0905 


0.0891 


0.0889 


0.2 


4.41 


13.5 


0.1703 


0.2947 


0.2872 


0.2866 


0.2 


4.48 


11.5 


0.2828 


0.4952 


0.4809 


0.4800 


0.2 


4.55 


9.9 


0.3877 


0.6558 


0.6371 


0.6363 


0.1 


4.61 


8.5 


0.4817 


0.7727 


0.7529 


0.7523 


0.1 


4.66 


7.3 


0.5641 


0.8534 


0.8348 


0.8344 


0. 1 


4.70 


6.3 


0.6352 


0.9071 


0.8910 


0.8907 


0.0 


4.74 


5.4 


0.6959 


0.9419 


0.9288 


0.9286 


0.0 





C 




25, o - eV «* = 


2.50 










Po 


P 


P 


Pc . 




1.77 


41.0 


0.0573 


0 . 0660 


0.0652 


0.0648 


0.5 


1.82 


37.1 


0.2137 


0.2733 


0.2663 


0.2638 


0.9 


1.86 


34.1 


0.3706 


0.4822 


0.4674 


0.4631 


0.9 


1.90 


31.7 


0.5076 


0.6497 


0.6299 


0.6249 


0.8 


1.93 


29. 7 


0.6206 


0.7710 


0.7497 


0.7451 


0.6 


1.95 


28.0 


0.7107 


0.8539 


0.8339 


0.8300 


0.5 


1.98 


26.5 


0.7812 


0.9085 


0.8912 


0.8881 


0.3 


1.99 


25.3 


0.8356 


0.9435 


0.9295 


0.9272 


0.2 
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optimal dispersion is not so different from the true 

maximum. The relative loss is at most 1%, and is 

2 2 

very small when a^/a is large. 

2 

The top curve in Figure 5 is the relative error Ao^ of 

the classical approximate formula, where the abscissa should 

be read as in this case. In spite of its simple form, it 

is indeed an excellent approximation, in particular for large 

2 

z 's. The z -dependence of the relative error Aa ' is, how- 
c c c 

ever, more sharp than ours. The middle curve in Figure 5 will 
be discussed in the next section. 

B. IMPROVEMENT OF THE APPROXIMATION 

In the preceding section, we investigated the accuracy 
of the approximate formula given by (4.10), and it seemed 
that (4.10) gave a lower bound to the true optimal ballistic 
dispersion. The reason it gives a possible lower bound, and 
a method of obtaining improved approximations will be con- 
sidered now. 

The Walsh optimal f 2 and others in (3.16) are functions 
of X, 

X = C(o^ + a)2^aV2 , (4.15) 

for the Simple Salvo MOdel, and so let us denote them as 

f 2 (x) instead of ' etc. With this notation, (3.18) 

is 
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p*(x)/D(x) . 



(4.16) 



2 2 

Equation (3.27) by which the approximate a /a is determined 
reads 





Pq (x)/Dq (x) 



(4.17) 



e 



and 



D(x) 




(4.18) 



because the damage function is assumed to be Gaussian. 

The function p*(x) is the double Fourier transform of 
the Walsh optimal p function given in Corollary 3.2 with 



given by (3.26) . Figure 6 illustrates what the curve of p* (x) 
is like. The curve decreases as ic is increased from zero, 
changes its sign, and then swings back to positive with small 
amplitude, and so forth. On the other hand, the approxi- 
mating function Pq C x) is equal to p* (x) at x = 0, has the 
same tangent there, and decreases exponentially. Pq (x) is 
larger than p*(x) for a wide interval of x as is seen in 
Figure 6. 



and Pq (-x) is an approximate formula for it 



Pq(x) > p* (x) 



0 < X < X 



0 * 
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p (x )/2Tra 




X 

Figure 6'. 
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Therefore, from (4.16) and (4.17), we have 



2 , 2 

-xa /a * 

e ^ > f^(x) , 0 < X < Xq 

2 

which suggests that is smaller than the true optimal 
value . 

One of the simplest ways to get an approximate formula 

2 

would be to define a, 

b 

2 2 
-0xa^/a 

e ^ = Pq(x)/Dq(x) (4.19) 

with 0 < 9 < 1 , or 

% " ' 0 < 6 < 1 , (4.20) 

and to determine the parameter 9 empirically to achieve 

2 

satisfactory agreement between the approximate and the 

2 2 
true 0 2 > Note that the definition of corresponds to 

approximating the p*(x) function by p^^(x) with an exponent 

1/9 times the exponent of Pq (x) . 

For example, let 1/9 be 1.15. Then, 

= 1.15 . (4.21) 

In Table 9A, Table 9B, and Table 9C, we present the accuracy 
investigation of (4.21). Accuracy is found much improved. 
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n 

5 

10 

15 

20 

25 

30 

35 

40 

n 

4 

8 

12 

16 

20 

24 

28 

32 

n 

2 

4 

6 

8 

10 

12 

14 

16 



Table 9 Accuracy of Approximation (4.2/) 



A z = 400, 9.89 







Po 


P 


P 


Pb 


^?‘/c 


9.34 


5.9 


0.0225 


0.0414 


0.0410 


0.0410 


0.0 


9.46 


4.5 


0.0634 


0.1484 


0.1456 


0.1455 


0.0 


9.57 


3.3 


0.1088 


0.2777 


0.2709 


0.2709 


0.0 


9.67 


2.3 


0.1556 


0.4063 


0.3952 


0.3952 


0.0 


9.76 


1.3 


0.2024 


0.5229 


0.5081 


0.5031 


0.0 


9.84 


0.5 


0.2484 


0.6230 


0.6056 


0.6056 


0.0 


9.92 


-0.3 


0.2931 


0.7061 


0.6873 


0.6873 


0.0 


9.98 


-1.0 


0.3362 


0.7733 


0.7542 


0.7542 


0.0 







B z = 100 


. o~b = 


4.47 










P« 


P 


? 


Pb 


^P/o 


4.31 


3.6 


0 . 0600 


0.0905 


0.0891 


0.0691 


0.0 


4.41 


1.4 


0.1703 


0.2947 


0.2872 


0.2872 


0.0 


4.48 


-0.4 


0.2828 


0.4952 


0.4809 


0.4809 


0.0 


4.55 


-1.8 


0.3877 


0.6558 


0.6371 


0.6371 


0.0 


4.61 


-3.1 


0.4817 


0.7727 


0.7529 


0.7528 


0.0 


4.66 


-4.1 


0.5641 


0.8534 


0.8348 


0.8347 


0.0 


4.70 


-5.0 


0.6352 


0.9071 


0.8910 


0.8908 


0.0 


4.74 


-5.8 


0.6959 


0.9419 


0.9288 


0.9286 


0.0 







o 

N 

It 


25, = 


1.76 








4 


Po 


P 


P 


Pb 




1.77 


-0.6 


0.0573 


0.0660 


0.0652 


0.0652 


0.0 


1.82 


-3.3 


0.2131 


0.2733 


0.2663 


0.2663 


0.0 


1.86 


-5.4 


0.3706 


0.4822 


0.4674 


0.4673 


0.0 


1.90 


-7.1 


0.5076 


0.6497 


0.6299 


0.6296 


0. 1 


1.93 


-8.6 


0.6206 


0.7710 


0.7497 


0.7493 


0.1 


1.95 


-9.7 


0.7107 


0.8539 


0.8339 


0.8333 


0. 1 


1.98 


-10.8 


0.7812 


0.9085 


0.8912 


0.8906 


0.1 


1.99 


-11.6 


0.8356 


0.9435 


0.9295 


0.9239 


0.1 
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2 

The maximum absolute relative error Ao, reduces to 12%, and 

b 

the relative loss in the salvo kill probability is, at most 
only 0.07%. The middle curve in Figure 5 is the relative 

A 2 

error Ao. . 
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V. SALVO MODEL WITH NON GAUSSIAN DAMAGE FUNCTION 



In this chapter we deal with a salvo kill model with 
circular normal errors and circularly symmetric general 
damage function. Under this general assumption, calculation 
of the salvo kill probability is very difficult, and time- 
consuming computation is needed for it. Here we propose an 
approximate method for calculating the salvo kill proba- 
bility. By this method computation is very easy and yet 
fairly high accuracy is expected. 

Assumptions in this chapter are as follows: A salvo of 

n weapons is fired. The impact point error is composed of 
two parts, the random bias common to all n weapons and the 
roxind-to-round error. We adopt a Cartesian coordinate sys- 
tem (x,y) the origin of which coincides with the center of 
impact of the n weapons. The random position of the target 
is denoted as (U,V) with respect to this coordinate system, 
and (U,V) is assumed to be circular normal centered at C0,0) 
with variance a^: 



f- Cu,v) = — e ^ . (5.1) 

^ 2tto^ 

The impact points of the ith weapon (X^,Y^) are independent 
and identically distributed circular normal random variables 
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2 

centered at (0,0) with variance q 



£2 (x,y) 





(5.2) 



The conditional kill probability of the target by a single 
weapon given that the target is at (u,v) and the weapon impact 
point is (x,y) is called as the damage function D. In this 
chapter, D is assumed to be a function of the miss distance 
r = >/(u-x) ^ + (v-y) ^ only. 



D = D(r) 



2 

r 



(u-x) ^ + (v-y) ^ 



(5.3) 



Further, the lethal area is denoted as A, and it is assumed 
that there exist finite "moments" of all degree as defined 
in (5.5), 



00 

2tt / D(r) rdr = A , (5.4) 

0 



r^ D(r) rdr < <» , j = 1,2,... (5.5) 

0 



The conditional kill probability by a single weapon 
given (U,V) being (u,v) is given by Cl. II, repeated here as 



p(u,v) = / / D(u-x,v-y) f 2 (x,y) dx dy , C5.6) 



E(r^) = 



2it 

A 
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and the salvo kill probability (1.2) is renamed as (5.7) . 

00 oo ^ 

P = 1 - / / {l-p(u,v)} f^(u,v) du dv . (5.7) 

— OO —00 

In the Simple Salvo Model in Chapter II with Gaussian 
assumption on the damage function, we have a simple expres- 
sion for p(u, v) ; 



2 2{a^+ah 

p(u,v) = ^ e ^ . (2.6) 

In the general model here, p(u,v) given by (5.6) is a com- 

r~2 2 

plicated function in vu +v , and computation of the kill 
probability is difficult. To overcome this difficulty, 
we propose an approximate formula for p(u,v), and then 
present a recursive method for computing the salvo kill 
probability with this approximate p(u,v) fianction. 

A. APPROXIMATE FORMULA FOR p(u,v) 

In this section, an approximate formula for the condi- 
tional kill probability p(u,v) is presented. As a corollary, 
we will get an approximation to the so-called Circular 
Coverage Function. 

Define the double Fourier transform of p(u,v) as 
- oo 00 i(o)^u+a)2v) 

pCw^,^^) = / / ® p(u,v) du dv , C5.8) 
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with similar notations being used for D(x,y) and f(x,y). 
Then, from (5.6), we have 

p(u)j^,W2) = D (ojj^ /0)2) f2 * (5.9) 

The Fourier transform of the circular normal density (5.2) 
is 



f 2 (o)^ , W 2 ) 





) 



(5.10) 



In the Fourier transform of the damage function, let the 
exponential function be expanded and integrated term by term. 



D (to ,( 1 ) 2 ) 



oo 00 i (to-, x+co^y) 
f f e D(x,y) dxdy 



00 00 CX3 



~ I f / '^ i^ (to, X + to^y) (x,y) dxdy 

j =0 -00 -00 J • 



By symmetry, the odd terms are zero, so 



k 00 00 



r \ ^ 1 21c 

D(to^,t02) = A I jjk TT A ^ (c0j_x+t02y) D(x,y) dxdy 

k =:0 * —00 —00 



^ ^ — 1 21c 21c 

= A y T / (/ (( 0 , cos 9 +to. sin 9 ) ^^de}r^^D(rIr dr 

k=0 '^00 ^ ^ 
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Therefore 



D(w, ,o)_) /A = 1 + I 
^ k=l 



(- 1 )^ 

uwr 



(2k-l) (2k-3) » ♦ . ^2k 
2k(2k-2) • • -2 



^ t 2^ 2> k 

) (w^+w^) 



(5.11) 



since 



(2k-l) (2k-3) ♦ ♦ .1 . > 



1 2k 

^ j (o)^ cos 9+^2 ~ 



2k(2k-2) • • -2 



, k i 1 



k = 0 



The first few terms in (5.11) are 



D(w^,o)2)/A = 



, 1 „ , 2, - 2 ^ 2. ^ 1 4- . 2 2. 2 

1 - j E(r ) (Uj^ ■'■“2^ ^ ^“1 "^“2^ 



Now let us define 



^a^‘^1'^2^ = A{1 



b, 2 , 2> . -c(wJ+W2) 

j(Wj^ + ui^) } e 



(5.12) 



and determine two constants b and c such that (5.11) and 
(5.12) coincide up to the fourth power of Wj^'s. Then we 
get 



b = I 7{E(r^) - E(r^)/2 , 



(5.13) 



and 
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c = -b/2 + E(r^)/4 . 

If (5.10) and (5.13) are siibstituted in (5.9), the result is 






= A{1 - 



b, 2 

2 ( 0^1 



s , 2^ 2, 
. 2uo~ 1 2 ’ 

+ W^) > e 



(5.14) 



where 



= O 2 - b + E(r^)/2 



(5.15) 



with b given by (5.13) . 

Next, the inverse transform of p (oj, ,U3_) is calculated: 

^ J. ^ 



p (u,v) 

d 



00 00 -i(uto,+vco_) 

■^/ / e 



(27T) 



— 00 —00 



Pa (^ 1 ,0^2) dWj^dw^ 



(5.16) 



The inverse transform of the first term in (5.14) 



Ae 



2,2 2 ^ 

-S (w^+w^) /2 



is easily calculated as 



2 , 2 

u +v 



27tS 



A ^ 2S' 



The second term in (5.14) is 
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Ab 2 2 " 2 

— (03^+002) e 



= -b 



d(S^) 



S / 2^ 2, 

+<^o) 

■{Ae ^ ^ ^ , 



and its inverse transform is given by 



2^ 2 
u +v 

V- d r A „ 2S^ , 

b ^ \ « e j 

d(S ) 2 ttS 



Therefore we have 



(u,v) 



2itS^ 



-U - -|-(1 - }e 

S 2S 



2S^ 



(5.17) 



It is noted that (5.17) is exact when the damage function 

is Gaussian with parameter a; (5.17) coincides with the exact 

2 

formula (2.6) because in this case A = 2ira , b = 0, and 




The cookie cutter damage function differs substantially 
from the Gaussian, so consideration of the former provides 
a good test of the robustness of (5.17) . The cookie cutter 
damage function is 



D(r) 




if 0 < r i a 



otherwise 



C5.18) 
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The probability 



p(u,v) 



00 00 

/ / D (u-x, v-y) f 2 (x,y) dx dy 

»-oo — oo 



2-no. 



(x-u) (y-v) ^ 

o 2 

2a_ 

j I e dx dy 

2 ^ 2 < 2 
X +y =a 



is called the Circular Coverage Function. In the literature, 

[~2 2 

the notation {s,/o .^,t:/o is used for it, where r = vu +v . 

A standard method of calculating the Circular Coverage Func- 
tion is the algorithm given by Brennan and Reed [Ref. 10] 
based on an infinite series expansion, but if the accuracy 

required is not too high, it is convenient to use closed form 

approximations. Standard handbooks [Refs. 11,12] recommend 

the next formula for ^ 



Pa(^/o2'^/c^2) 




(5.19) 



with 



v2 2 2,. 

S = 02 + a /4 



(5.20) 



According to (5.17), the Circular Coverage Function 

2 2 
(x-r) ^+y^ 

o 2 

1 ^^2 
P(a/o^,r/a^) = J J ® 



C5.21) 



2 2 2 . 2 

X +y la 
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can be approximated by 



P^(a/a2,r/a2) 




2S S 2S 



(5.22) 



where 



b 



a^/4/3 , = 02 + (1 - l//3)/4 



(5.23) 




cookie cutter damage function. 

In Table 10, the Circular Coverage Functio P(R,d) is 

tabulated for R = 0.1 (0.1) 2.0 and d=0.0 (0.5) 2.5. In 

Table 11, the error A^(R,d) = P^(R,d) - P(R,d) is shown. 

For small R, the error is very small; for instance the error 

is less than 10 when R 1 1.0. But when R is increased 

beyond 1, the error grows rapidly. 

In Table 12, the error A (R,d) = P (R,d). - P(R,d) of 

^ 0 . 

our approximate formula is given. The accuracy is found much 

better than the above-mentioned formula. The absolute error 

-4 < -2 

for d = 0 is less than 5 x 10 when R = 1.0, less than 10 

even when R = 1.6 and usually shows much higher accuracy. 

B. RECURRENCE FORMULA FOR THE SALVO KILL PROBABILITY 

If the conditional kill probability p(u,v) by a single 
weapon is given in the form (5.17), we can easily compute the 
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Table 10 



Circular Coverage Function P(R,d) 



d 



R 


0 . 00 


0.50 


1.00 


1.50 


2.00 


2.50 


0.1 


0.0050 


0.0044 


0.0030 


0.0016 


0.0007 


0.0002 


0.2 


0.0198 


0.0175 


0.0121 


0.0065 


0.0027 


0 . 0009 


0.3 


0.0440 


0.0389 


0.0270 


0.0146 


0.0062 


0.0021 


0.4 


0.0769 


0.0682 


0.0476 


0.0261 


0.0112 


0.0038 


0.5 


0.1175 


0.1045 


0.0735 


0.0408 


0.0179 


0.0062 


0.6 


0.1647 


0.1470 


0.1043 


0.0589 


0.0264 


0.0094 


0.7 


0.2173 


0.1946 


0.1397 


0.0803 


0.0369 


0.0135 


0.8 


0.2739 


0.2463 


0.1790 


0.1050 


0.0495 


0.0187 


0.9 


0.3330 


0.3009 


0.2217 


0.1328 


0.0645 


0.0252 


1.0 


0.3935 


0.3573 


0.2671 


0.1638 


0.0819 


0.0332 


1.1 


0.4539 


0.4144 


0.3146 


0.1976 


0.1019 


0.0428 


1.2 


0.5132 


0.4712 


0.3635 


0.2341 


0.1247 


0.0543 


1.3 


0.5704 


0.5267 


0.4132 


0.2730 


0.1501 


0.0679 


1.4 


0.6247 


0.5802 


0.4628 


0.3138 


0.1784 


0.0838 


1.5 


0.6753 


0.6309 


0.5120 


0.3563 


0.2092 


0.1021 


1.6 


0.7220 


0.6785 


0.5599 


0.3999 


0.2426 


0.1228 


1.7 


0.7643 


0.7224 


0 . 6062 


0.4442 


0.2784 


0.1462 


1.8 


0.8021 


0.7625 


0.6504 


0.4887 


0.3161 


0.1723 


1.9 


0.8355 


0.7987 


0.6921 


0.5329 


0.3556 


0.2010 


2.0 


0.8647 


0.8309 


0.7310 


0.5763 


0.3965 


0.2321 
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Table 11 Error Ap ^ ( R *< 3 ) 



d 



R 


0 . 00 


0.50 


1.00 


1.50 


2.00 


2.50 


0.1 


— 0 . 0000 


— 0 , 0000 


— 0 . 0000 


— 0 , 0000 


— 0 . 0000 


— 0 . 0000 


0.2 


0 . 0000 ' 


0 . 0000 


0 • 0000 


— 0 , 0000 


— 0 . 0000 


- 0.0000 


0.3 


0 . 0000 


0 . 0000 


0 • 0000 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


0.4 


0.0000 


0 . 0000 


0 • 0000 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


0.5 


0.0001 


0.0001 


0 • 0000 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


0.6 


0 . 0004 


0 . 0003 


0 • 0000 


- 0.0001 


— 0 . 0001 


— 0 . 0000 


0.7 


0.0010 


0 . 0007 


0 • 0001 


- 0.0002 


- 0 . 0002 


— 0 . 0000 


0.8 


0.0020 


0.0014 


0 • 0003 


— 0 . 0004 


— 0 . 0003 


— 0 . 0001 


0.9 


0.0038 


0.0027 


0 • 0005 


- 0.0007 


- 0.0006 


— 0 . 0002 


1.0 


0.0065 


0 . 0046 


0.0010 


- 0.0012 


- 0.0011 


— 0 . 0003 


1.1 


0.0106 


0.0076 


0.0018 


- 0.0018 


- 0.0019 


- 0.0007 


1.2 


0.0162 


0.0118 


0 . 0030 


- 0.0026 


— 0 . 0030 


- 0.0011 


1.3 


0.0236 


0.0174 


0 • 0048 


- 0.0036 


- 0.0045 


- 0.0019 


1.4 


0.0330 


0.0246 


0.0074 


- 0.0047 


- 0.0065 


- 0.0030 


1.5 


0.0447 


0.0337 


0.0109 


- 0.0058 


— 0 . 0090 


— 0 . 0046 


1.6 


0.0585 


0.0447 


0.0155 


- 0.0068 


- 0.0121 


- 0.0067 


1.7 


0.0746 


0.0578 


0.0213 


- 0.0076 


- 0.0157 


- 0.0095 


1.8 


0.0929 


0.0728 


0.0286 


— 0 . 0080 


- 0.0197 


- 0.0131 


1.9 


0.1132 


0.0898 


0.0374 


- 0.0077 


- 0.0240 


- 0.0174 


2.0 


0.1353 


0.1086 


0.0478 


- 0.0065 


- 0.0286 


- 0.0225 



81 



Table 12 



Error ( R , d ) 



d 



R 


0.00 


0.50 


1.00 


1.50 


2.00 


2.50 


0.1 


— 0 , 0000 


- 0 . 0000 


- 0.0000 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


0.2 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


- 0.0000 


- 0 . 0000 


0.3 


- 0.0000 


- 0.0000 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


— 0 . 0000 


0.4 


- 0.0000 


— 0 . 0000 


— 0 . 0000 


0 . 0000 


— 0 . 0000 


- 0.0000 


0.5 


— 0 . 0000 


— 0 . 0000 


0 • 0000 


0 . 0000 


0 . 0000 


- 0 . 0000 


0.6 


- 0 . 0000 


— 0 . 0000 


0 • 0000 


0 . 0000 


0 . 0000 


- 0 . 0000 


0.7 


— 0 . 0000 


- 0 . 0000 


0 • 0000 


0 . 0000 


0.0000 


— 0 . 0000 


0.8 


- 0.0001 


- 0.0000 


0 • 0000 


0 . 0000 


0 . 0000 


- 0.0000 


0.9 


- 0.0001 


— 0 . 0001 


0 . 0000 


0 . 0000 


0 . 0000 


— 0 . 0000 


1.0 


- 0 . 0003 


- 0.0002 


0 • 0000 


0.0001 


0 . 0000 


- 0 • 0000 


1.1 


- 0.0006 


— 0 . 0003 


0 . 0000 


0 . 0001 


0 • 0000 


- 0 . 0000 


1.2 


- 0.0011 


- 0.0007 


0.0000 


0 . 0003 


0 . 0001 


- 0 . 0000 


1.3 


- 0.0019 


- 0.0012 


0 • 0000 


0 . 0005 


0.0002 


- 0.0001 


1.4 


- 0.0031 


- 0.0019 


0 . 0000 


0.0008 


0.0003 


- 0.0001 


1.5 


— 0 . 0049 


- 0.0031 


- 0 . 0000 


0.0012 


0 . 0006 


— 0 . 0001 


1.6 


- 0.0075 


- 0.0048 


- 0 . 0001 


0.0019 


0 • 0009 


- 0 . 0002 


1.7 


- 0.0110 


- 0.0072 


— 0 . 0003 


0.0027 


0.0015 


- 0.0002 


1.8 


- 0.0157 


- 0.0104 


— 0 . 0007 


0.0039 


0.0023 


- 0.0002 


1.9 


- 0.0218 


- 0.0146 


- 0.0012 


0.0053 


0.0034 


— 0 . 0001 


2.0 


- 0.0294 


- 0.0199 


- 0.0021 


0.0071 


0.0049 


0 . 0001 
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salvo kill probability in a recursive way. We first observe 
that the salvo kill probability 



OO 00 

P = 1 - / / {1 - p (u, v) }^f, (u, v) du dv 

— OO — OO 



where f^(u,v) is the circular normal density (5.1) and 

2 

r 

2 ~ 2 

p(u,v) = (y+3 e , r^= u^+v^ , (5.24) 

2S^ 

can be expressed as 

P = 1 - K(n,0) , (5.25) 



where 

OO 

K(n,0) = / {1 - (y + 3t/p)e"'‘^/'^}^e"^ dt , (5.26) 

0 

2 2 

and where p = S /o^. This is simply a matter of substitution 
and the introduction of polar coordinates. We then have 

Proposition 5.1 
The function 



K(n,j) = / {1 - (y + 3t/p)e"^/P}^e ^dt (5.27) 

0 
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which has two integer arguments and three parameters y, 6, 
and p satisfies the following recurrence relation; 



K(n,j) = 


n=l,2, ... 7 j=0,l,... (5.28) 


with 





K(0,j) = / j = 0,1, ... (5.29) 



Proof 




From 


(5.27) , 


K(n,j) = 


/“{I - it 

0 

0 



The first term in the right-hand-side is K(n-l,j). The second 
terms is, when integrated by parts. 
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00 



-/ {1 - (u + 6t/p) 4 6t/p) 

0 



= - ^{1- (y+6t/p)e-^/P}^e-(3/P-^l)^f 

0 



- ^(j/p+1) / a - (U + 0t/p)e"^/P}^e" ^ dt 

0 

- sfu - (y + 6t/p)e-'^/P}^e-^/Pe-(^/P-^l^^ dt 

0 

= ^(1-y)"^ - i^K(n,j) - 6K(n-l,j4-l) . 

Therefore 

K(n,j) = K(n-l,j) + ^(1-y)" - K(n,j) - 3K(n-l,j4-l) 

from which (5.28) is immediate. From the definition of 
K(0,j), it is clear that (5.29) holds. 

Note 1 

We use the notations y and p here in the same way as 
in Chapter II without any subscripts because no confusion 
is expected. Their definitions are different, but they are 
just the counterparts of y and p defined in Chapter II. 

Note 2 

It is worthwhile to note that the recurrence relation 
(5.28) has a similar form to (2.17) . The p function (5,24) 
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involves an extra parameter 6 which was missing in the Simple 
Salvo Model in Chapter II, and therefore, it is necessary 
to introduce another index j . 
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VI, COOKIE CUTTER SALVO MODEL 



Consider that n weapons are fired at a target in a 

salvo. The impact point error is composed of two parts, 

the random bais and the dispersion. We adopt a Cartesian 

coordinate system (x,y) the origin of which coincides with 

the center of impact of the n weapons. The random position 

of the target is denoted as (U,V) with respect to this 

coordinate system, and (U,V) is assumed to be circular normal 

2 

centered at (0,0) with variance 



(u,v) 



2iTaJ 



2 ^ 2 
u +v 



2a: 



( 6 . 1 ) 



The impact points of the ith weapon (X^,Y^) are indepen- 
dent and identically distributed circular normal random 

2 

variables centered at (0,0) with variance 02 • 



f2 (^/Y) 





( 6 . 2 ) 



The kill probability by a single weapon conditional on 
(U,V) = (u,v) and (X,Y) = (x,y) is a function of u-x and 
v-y, D(u-x,v-y) , and in this chapter it is assxamed that 
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D (u-x, v-y) 



1 if (u-x)^ + (v-y) ^ 
0 otherwise , 



with lethal area 



(6.3) 



^ 2 
A = TT a . 

In the following, we call this model the Cookie Cutter 
Salvo Model. 

The conditional probability that the target is destroyed 
by a single weapon given the random position of the target 
being (u,v) is 



00 oo 

p(u,v) = / / D(u-x, v-y) f 2 (x,y) dx dy , 

— oo —00 



and with (6.2) and (6.3) , we have 



p(u,v) = 



/ / 



(u-x) (v-y) ^ 

202 

e dxdy = P(a/a, 



,v^u^+v^ /a_) , 



x^+y2<a2 



(6.4) 



where P(R,d) is the Circular Coverage Function, 
The salvo kill probability 



OO 00 



P = 1 - / / {1 - p(u,v) }^f j^(u, v) du dv , 



—00 — oo 
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with (6.1) and (6.4) is readily rewritten as 



oo 



p 



c 



1 - / {1 - P (a/02 ^ /^cj^/a^) }"e~^ dt . 



0 



(6.5) 



We carried out a number of computations of (6.5) with 



various sets of parameters a^r a, and n by the Simpson 



method. A robust algorithm for computing the Circular 
Coverage Function P(R,d) for wide ranges of R and d was needed. 
For this purpose, we used the Brennan-Reed formula in a modi- 
fied form. 

Brennan-Reed Formula: 



OO 



P(R,d) = I g k 
n=0 ^ " 




d^/(2n) / n 



1,2 




e 




g. 



n 




n 



1,2 
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Modified Formula: 



P(R,d) = e 



2 

-d /2 



I 

n=0 



n n 



(6.7) 



Wl = dV(2n) , n = 1,2, .. . 



h 



0 



1 , 



C 



n 



=11-1 * ^n' 



n 



1 , 2 , 






( 6 . 8 ) 



a /a , 
n'^ n-1 



= RV(2n) 



n = 1,2, 



The modified formula is especially of use in preventing 
underflow possible in cases with large R and d. 

In their review paper, Eckler and Burr presented a short 
table of the salvo kill probability of the Cookie Cutter 
Salvo Model. But the accuracy is, as they admitted, sup- 
posedly low because it is obtained by Monte Carlo simula- 
tion. In Table 13, we give the correct values. The figures 
in parentheses are from the table by Eckler and Burr. 

In the following sections, comparison of the Cookie Cutter 
Salvo Model and the Simple Salvo Model is dealt with and 
then the optimal ballistic dispersion is investigated, where 
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Table 13 Salvo Kill Probability ; Kisi/( Eckler-Burr ) 



/ cr. 





0.25 


0.50 


0.75 


0.6 


0.221 


0.228 


0.199 




(0.224) 


(0.235) 


(0.201 ) 


1.2 


0.573 


0.587 


0.557 




(0.574) 


(0.582) 


(0.563) 


1.8 


0.834 


0.836 


0.813 




(0.834) 


(0.832) 


(0.814) 


2.4 


0.954 


0.952 


0.939 




(0.949) 


(0.951) 


(0.942) 


3.0 


0.991 


0.990 


0.985 




(0.990) 


(0.988) 


(0.985) 



B n = 5 



C n = 20 



Ox/cr^ 



Q /<7, 


0.50 


0.75 


1.00 


0.4 


0.231 


0.213 


0.177 




(0.225) 


(0.214) 


(0.187) 


0.8 


0.525 


0.551 


0.514 




(0.521) 


(0.551) 


(0.516) 


1.2 


0.736 


0.772 


0.764 




(0.728) 


(0.773 ) 


(0.763) 


1.6 


0.872 


0.894 


0.896 




(0.865) 


(0.892) 


(0.893 ) 


2.0 


0.946 


0.957 


0.958 




(0.945) 


(0.956 ) 


(0.957) 







0“ i/ <Tf 






<3/0-, 


0.50 


0.75 


1.00 


1.50 


0.2 


0.221 


0.210 


0.176 


0.115 




(0.209) 


(0.201) 


(0.176) 


(0.113) 


0.4 


0.458 


0.525 


0.504 


0.381 




(0.447) 


(0.531 ) 


(0.499) 


(0.386) 


0.6 


0.605 


0.713 


0.740 


0.648 




(0.595) 


(0.716) 


(0.740) 


(0.640) 


0.8 


0.715 


0.818 


0.861 


0.828 




(0.706) 


(0.817) 


(0.846) 


(0.828) 


1.0 


0.801 


0.884 


0.922 


0.923 




(0.792) 


(0.882) 


(0.924) 


(0.924) 


1.2 


0.865 


0.928 


0.956 


0.966 




(0.856) 


(0.926) 


(0.958) 


(0.970) 


1 .4 


0.912 


0.956 


0.975 


0.985 




(0.900) 


(0.953) 


(0.975) 


(0.984) 
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the approximate formula presented in Chapter III and the 
approximate method of computing the salvo kill probability 
developed in Chapter V are found to be useful. 

A. COOKIE CUTTER VS. SIMPLE SALVO MODEL 

The Cookie Cutter Salvo Model differs from the Simple 
Salvo Model of Chapter II only in the damage function. In 
Chapter II, we assumed 

(u-x) (v-y) ^ 

D( u-x, v-y) = e ^ , 

2 

with lethal area 2-na . For comparison's sake, the lethal 
area of the two salvo models are now set equal 

= 2a^ . (6.9) 

The two models share the common parameters a, and n, 

only being different with respect to the shape of the damage 
function. The quantities associated with the Cookie Cutter 
Salvo Model will be given suffix c in the following when 
necessary. 

It is interesting to compare the salvo kill probabilities 
of the two models. We begin with three extreme cases in which 
the salvo kill probability has closed form expressions. 

(1) The case n = 1. 

The formula (6.5) with (6.4) is easily integrated and 

give 
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( 6 . 10 ) 



-a^/(ol+oh 

P = 1 - e -L ^ 

c 

The corresponding single shot kill probability under the 
assumption of the Gaussian damage function is 

P = 1 ^ ^ ^ . (6.11) 

l+aV(a£ + a‘) 

Proposition 6.1 

The Cookie Cutter Salvo Model with n = 1 gives a higher 
single shot kill probability than that of the Simple Salvo 
Model . 



P > P . (6.12) 

c 

Proof 

Recall that e ^ < l/(l+x) for x > 0. Let a^/Ca^+a^) = x 
in (6.10) and (6.11), then (6.12) is obvious. □ 

Generally, in the salvo with n > 1, this inequality 
does not necessarily hold as will be shown in the following. 

(2) The case = 0 

Let ^ 0 in (6.5) , then we have 

2 , 2 
-na /a, 

P = 1 - e . (6.131 

c 
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The corresponding formula in the Simple Salvo Model is given 
by (2 . 23) , namely 



P = 1 - 1/(1 + 0 ^/ 0 ^)'^ . 



(6.14) 



(3) The case ~ ^ 

Let 0 in (6.4) . We get 



P 



c 



1 



-a 



e 




(6.15) 



The corresponding formula in the Simple Salvo Model is given 
by (2.24) : 



P 



n (n— 1) • • • 1 

(n+A) (n-l+A) • • • (1+A) ' 



(6.16) 



with 



2 , 2 
= a /o^ 



Proposition 6.2 

Between the salvo kill probabilities of the Cookie 
Cutter Salvo Model and the Simple Salvo Model with the same 
parameters, the following relationships hold; 

If = 0 , then 



P > P . 
c 



(6.17) 
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If O 2 - 0, then 

< P iff o^/a^ < 1/Aq , (6.18) 

where Aq is a positive solution to the equation 



n 

A - ^ ln(l + A/j) = 0 . (6.19) 

j=l 



Proof 

The proof of (6.17) is just the same as that of Proposi- 
tion 6.1. As to the second assertion, let 

P = 1 - 

in (6.16). Then 

n n 

g(A) = I ln(l + A/j) , dg/dA = I l/(j+A) . 

j=l j=l 

2 2 

since g(0) =0, dg/dA | > 1, d g/dA < 0, we have 

A=0 

g(A) 5 A iff A $ Aq . 

where A^ is given by (6.19). The relation (6.18) is 
immediate. □ 
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As to the general case, we must compare the two kill 

probabilities numerically. To illustrate the general idea. 

Table 14 is given. In Table 14A, for n = 4 is tabulated, 

and the corresponding Simple Salvo Model kill probability 

P is in Table 14B. It is observed that P 's in the column 

c 

2 2 

Oj^/a = 0 are larger than P's as stated in Proposition 6.2. 

2 2 

P^ ' s in the first row corresponding to <^ 2 ^^ = 0 is smaller 

2 2 

than P's except for the case Oj^/a = 0, since 1 /Aq = 0.216 
for n = 4 . 

As a general tendency, P is larger than P when is 
relatively larger than P^ is smaller than P when O 2 

is considerably smaller than as is seen in the upper 

right corner of the table. 

B. OPTIMAL BALLISTIC DISPERSION 

We carried out a golden search calculation to get the 
2 2 

optimal ratio <^ 2 ^^ ^ number of cases with different 

2 2 

combinations of n and o^/a. . The results are tabulated in 

2 2 

Table 15. Table 15A gives the optimal ^ 2^/01 of the Cookie 

Cutter Salvo Model, and Table 15B gives the associated 

salvo kill probability, namely the maximum salvo kill proba- 

2 2 

bility for a given set of n and a-^/a. . It is noted that this 
calculation is very time-consuming. The corresponding 
values for the Simple Salvo Model are tabulated in Table 16. 

Before starting comparison of the optimal ballistic 
dispersion of the Cookie Cutter Salvo Model with that of the 
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Table 14 



Salvo Kill Probability, n = 4 



A Cookie Cutter Salvo Model, P<. 



<r,V «. * 



aC/c(2 


0.0 


2.0 


4.0 


6.0 


• 

CD 


10.0 


0.0 


1.000 


0.393 


0.221 


0.154 


0.118 


0.095 


0.5 


1.000 


0.605 


0.389 


0.285 


0.224 


0.185 


1.0 


0.982 


0.626 


0.425 


0.320 


0.256 


0.214 


1.5 


0.931 


0.615 


0.435 


0.335 


0.272 


0.229 


2.0 


0.865 


0.592 


0.432 


0.339 


0.279 


0.237 


2.5 


0.798 


0.565 


0.423 


0.337 


0.281 


0.240 


3.0 


0.736 


0.538 


0.412 


0.333 


0.279 


0.241 


3.5 


0.681 


0.511 


0.399 


0.326 


0.276 


0.240 


4.0 


0.632 


0.486 


0.385 


0.319 


0.272 


0.238 







B 


Simple Salvo 


Model , 


P . 










cr,*/a 








/«* 


0.0 


2.0 


4.0 


6.0 


8.0 


10.0 


0.0 


1 .000 


0.594 


0.382 


0.280 


0.221 


0.183 


0.5 


0.988 


0.621 


0.421 


0.317 


0.253 


0.211 


1.0 


0.938 


0.613 


0.431 


0.331 


0.268 


0.225 


1.5 


0.870 


0.589 


0.427 


0.334 


0.274 


0.232 


2.0 


0.802 


0.560 


0.417 


0.331 


0.274 


0.233 


2.5 


0.740 


0.531 


0.404 


0.325 


0.271 


0.233 


3.0 


0.684 


0.503 


0.390 


0.317 


0.267 


0.230 


3.5 


0.634 


0.476 


0.375 


0.309 


0.262 


0.227 


4.0 


0.590 


0.452 


0.361 


0.300 


0.256 


0.223 
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Table 15 Optimal Ballistic Dispersion and 
Maximum Probability ( Cookie Cutter Salvo Model ) 



A Optimal Ratio d'ac* / at ^ 








2 


4 


8 


16 


32 


64 


2 


0.4 


1.0 


l.T 


2.7 


4.1 


6.1 


4 


0.8 


1.6 


2.6 


3.9 


5.8 


8.6 


6 


1.2 


2.1 


3.2 


4.9 


7.1 


10.5 


8 


1.6 


2.5 


3.8 


5.6 


8.3 


12.0 


10 


1.8 


2.9 


4.4 


6.4 


9.3 


13.4 


12 


2.1 


3.3 


4.9 


7.0 


10.1 


14.7 



B Maximum Salvo Kill Probability 

p 

C max 



n 





2 


4 


8 


16 


32 


64 


2 


0.476 


0.626 


0.782 


0.905 


0.973 


0.996 


4 


0.300 


0.435 


0.601 


0.769 


0.899 


0.971 


6 


0.224 


0.339 


0.492 


0.667 


0.826 


0.935 


8 


0.182 


0.281 


0.420 


0.590 


0.762 


0.897 


10 


0.154 


0.241 


0.368 


0.531 


0.708 


0.859 


12 


0.135 


0.213 


0.329 


0.484 


0.661 


0.823 
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Table 16 Optimal Ballistic Dispersion and 
Maxmum Probability ( Simple Salvo Model ) 



A Optimal Ratio (Tz* /o(^ . 



n 





2 


4 


8 


16 


'32 


64 


2 


0.5 


0.6 


1.3 


2.2 


3.6 


5.6 


4 


0.5 


1.1 


2.0 


3.4 


5.3 


8.1 


6 


0.5 


1.4 


2.6 


4.3 


6.6 


10.0 


8 


0.7 


1.8 


3.1 


5.0 


7.7 


11.5 


10 


0.9 


2.1 


3.6 


5.7 


8.7 


12.9 


12 


1.0 


2.3 


4.0 


6.3 


9.5 


14.1 



B Maximum Salvo Kill Probability 

P ma X 



n 





2 


4 


8 


16 


32 


64 


2 


0.450 


0.622 


0.780 


0.904 


0.973 


0.996 


4 


0.293 


0.431 


0.599 


0.768 


0.899 


0.971 


6 


0.217 


0.334 


0.489 


0.665 


0.826 


0.935 


8 


0.173 


0.274 


0.416 


0.588 


0.761 


0.896 


10 


0.144 


0.233 


0.363 


0.528 


0.707 


0.858 


12 


0.124 


0.204 


0.323 


0.481 


0.659- 


0.823 
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Simple Salvo Model, we recall the approximate formula (3.22) 
given in Chapter III. 



o oo 00 

/ / x^{p*(x,y) - D(x,y) }dx dy ,0]'*’ . (6.20) 

—00 — OO 



In the case of the Cookie Cutter Salvo Model, 



00 oo 

^ D(x,y) dx dy = a^/2 , 

— oo —00 



whereas in the Simple Salvo Model, 



oo oo 

^ / / x^ D(x,y) dx dy = . 

—00 —00 



Therefore, (6.20) states that the optimal ballistic dispersion 

of the Cookie Cutter Salvo Model is larger than that of the 

2 

Simple Salvo Model approximately by a /2: 

- a^/a^ = 0.5 . (6.21) 

Now we subtract entities in Table 16A from those in Table 

2 2 2 2 

15A and get the difference - o^/o. which is tabulated in 

Table 17A. The values are positive with only one exception, 
and the predicted value 0.5 is observed in the upper right 
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Table 17 Comparison of Cookie Cutter Salvo Model 

and Simple Salvo Model 



A Difference of Optimal Ratio <Tic/o(^- 



n 



cr,^/o(" 


2 


4 


8 


16 


32 


64 


2 


-0.1 


0.4 


0.4 


0.5 


0.5 


0.5 


4 


0.3 


0.5 


0.6 


0.5 


0.5 


0.5 


6 


0.7 


0.7 


0.6 


0.6 


0.5 


0.5 


8 


0.9 


0.7 


0.7 


0.6 


0.6 


0.5 


10 


0.9 


0.8 


0.8 


0.7 


0.6 


0.5 


12 


1.1 


1.0 


0.9 


0.7 


0.6 


0.6 



B ■ Difference of Maximum P's 

“ P max 



n 





2 


4 


8 


16 


32 


54 


2 


0 .026 


0.004 


0.002 


0.001 


0.000 


0 . 000 


4 


0.007 


0.004 


0.002 


0.001 


0.000 


0.000 


6 


0.007 


0.005 


0.003 


0.002 


0.000 


0.001 


8 


0 . 009 


0.007 


0.004 


0.002 


0.001 


0.001 


10 


0.010 


0.008 


0.005 


0.003 


0.001 


0.001 


12 


0.011 


0 . 009 


0 . 006 


0.003 


0.002 


0.001 
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corner of the table. In the below left corner of the table, 

the figures are not close to the value 0.5, but these cases 

correspond to low salvo kill probability. 

In Table 17B, the difference of the maximum probabilities 

of each model is tabulated. The figures are all positive, 

and lead to the speculation that P > P in all cases 

^ c max max 

despite the fact that P^ may be smaller than P when is 
not chosen optimally. The author has not yet succeeded in 
obtaining a proof of this speculation, but it might be the 
counterpart for n > 1 of Proposition 6.1 for n = 1. 

In Chapter V, we developed an approximate method of 
computing the salvo kill probability which can be applied 
to the Cookie Cutter Salvo Model. The salvo kill proba- 
bility is given by 



P 



1 - K(n,0) 



(5.25) 



where 



00 



"■t/p^n^-( j/p+1) ^ dt 



K(n,j) = / {1 - (y+3t/p)e 

0 



(5.27) 



satisfies the recurrence relation 




n+j+p 



•K(n-l,j+l), (5.28) 



n 



1,2 



j = 0,1 
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wi th 



K(0, j) 




(5.29) 



The parameters are given by (5.22), (5.23), and (5.24). 



(a^/S^) (1 - a^/2/3 S^) 



4 — 4 

B = a /2/3 S , 



( 6 . 22 ) 



S 



2 



+ a^(l - l//3)/2 . 



Using these formulae together with the golden section 
search calculation, we obtained Table 18. Table 18A gives 
the approximate optimal ballistic dispersion of the Cookie 
Cutter Salvo Model, and the associated salvo kill probability 
is in Table 18B. 

2 2 

In Table 19A, the difference between the optimal 

2 2 

given in Table 15A and the approximate Table 18A 

is tabulated. The difference decreases towards the upper 
right corner of the table. 

In Table 19 B, the salvo kill probability of the Cookie 
Cutter Salvo Model associated with the approximately optimal 
ballistic dispersion mentioned above is given. It corres- 
ponds to the maximum P given in Table 15B. In spite of the 

2 2 

discrepancies observed in Table 19A, the approximate o /a 

2 2 

gives almost the same kill probability as the true gives. 
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Table 18 Optimal Ballistic Dispersion 
and the Associated Salvo Kill Probability 
( Cookie Cutter Salvo Model, Approximation ) 



A Optimal Ratio • 



n 





2 


■ 4 


8 


16 


32 


64 


2 


0.5 


1 .0 


1.7 


2.7 


4.1 


6.1 


4 


0.7 


1.5 


2.5 


3.9 


5.8 


8.6 


6 


1.0 


1.9 


3.1 


4.8 


7.1 


10.5 


8 


1.2 


2.2 


3.6 


5.5 


8.2 


12.0 


10 


1.3 


2.5 


4.1 


6.2 


9.2 


13.4 


12 


1.5 


2.8 


4.5 


6.8 


10.0 


14.7 



B Maximum Salvo Kill Probability 



n 





2 


4 


8 


16 


32 


64 


2 


0.475 


0.625 


0.782 


0.905 


0.973 


0.996 


4 


0.297 


0.432 


0.600 


0.768 


0.899 


0.971 


6 


0.218 


0.334 


0.490 


0.665 


0.826 


0.935 


8 


0.173 


0.274 


0.416 


0.588 


0.762 


0.396 


10 


0.144 


0.234 


0.363 


0.528 


0.707 


0.358 


12 


0.124 


0.204 


0.323 


0.481 


0.659 


0.323 
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Table 19 Accuracy of the Approximation 
( Cookie Cutter Salvo Model ) 



A Error in the Optimal Ratio ^ 



n 



0*1 V c< ^ 


2 


4 


8 


16 


32 


64 


2 


-0.1 


0.0 


0.0 


0.0 


0.0 


0.0 


4 


0.1 


0.1 


0.1 


0.0 


0.0 


0.0 


6 


0.2 


0.2 


0.1 


0.1 


0.0 


0.0 


8 


0.4 


0.3 


0.2 


0.1 


0.1 


0.0 


10 


0.5 


0.4 


0.3 


0.2 


0.1 


0.0 


12 


0.6 


0.5 


0.4 


0.2 


0.1 


0.0 



B Salvo Kill Probability Associated with 
the Approximately Optimal Ballistic Dispersion. 

? C <Ta ) 



n 





2 


4 


8 


16 


32 


£4 


2 


0.475 


0.626 


0.782 


0.905 


0.973 


0.996 


4 


0.300 


0.435 


0.601 


0.769 


0.899 


0.971 


6 


0.223 


0.339 


0.492 


. 0.667 


0.826 


0.935 


8 


0.181 


0.280 


0.420 


0.590 


0.762 


0.397 


10 


0.153 


0.241 


0.368 


0.531 


0.708 


0.359 


12 


0.134 


0.212 


0.328 


0.484 


0.661 


0.323 
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VII. SALVO MODELS WITHOUT CIRCULAR SYMMETRY 



The approximate formula for the optimal ballistic dis- 
persion presented in Chapter III covers the salvo models 
without circular symmetry. In this chapter the formulae 
for computing the salvo kill probability of these models 
are dealt with. 

Suppose that a salvo of n weapons is fired against a 

target. The impact point error is composed of two parts, 

one being common to all the weapons, and another being 

round-to-round dispersion. The random target position (U,V) 

with respect to the coordinate system (x,y) is assumed to 

be elliptical normal centered at (0,0) with variances 
2 2 

a and a respectively. 



The impact points of the ith weapon (X^,Y^) are indepen- 
dent and identically distributed elliptical normal random 



u 



V 



f (u, v) 





( 7 . 1 ) 




f2 (x,y) 



1 




(. 7 . 2 ) 



2n O^Oy 
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The kill probability D by a single weapon conditional on the 
target being at (u,v) and the weapon impact point being (x,y) 
is a function of u-x and v-y only, and is given by the 



elliptical Gaussian function with parameters 
Specifically, 



a and a . 
X y 



(u-x) (v-y) ^ 

2 ~ 2 
2a 2a 

D(u-x,v-y) = e ^ ^ . (7.3) 

It might be an appropriate model, e.g., for a ground target 
vs. a weapon impacting the ground obliquely, because in this 
case the scattering of the fragments and thus the conditional 
kill probability D is by no means circular symmetric as 
as Slimed in the previous chapters. For brevity, we call 
this model the Elliptical Normal Salvo Model. 

The kill probability of a single weapon conditional on 
(U,V) being (u,v) is given by 



oo oo 



p(u,v) = / / D(u-x, v-y) f „ (x,y) dx dy , 



— OO — oo 



and using (7.2) and (7.3) we get 



u 



V 



p(u,v) = 






2 2 
2{a+a^) 

X X 



J, 2, 2. , 2^ 2, 



o / 2 . 2. 

2 (a +o ) 

y 



(7.4) 
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The salvo kill probability is then given by 



2 2 

1 °° °° 2a 2a 

P = 1 - 2 ir ' a o / / {l-p(u,v)}^e ^ dudv . (7.5) 

From this equation, Grxibbs derived the following formula: 
Proposition 7.1 (Grubbs) 

The kill probability of the Elliptical Normal Salvo 
Model stated above is given by 



P 



I (- 1 ) 

j=l ^ 



' PxPy 
(j+Pj,) (j + Py) 



(7.6) 



where 



y = a a / '/(a^ + a^) (a^ + a^) 
X y X X y y 



X 



/ 2 , 2 , , 2 

(a + a ) /a 

X X U 



(7.7) 



(a^ + a^) /a^ 
y Y ' V 



The Grubbs formula is not suited for calculation when n 
is large, since it is an alternating series. To overcome this 
difficulty, Breaux and Mohler [Ref. 13] gave a method for calcu- 
lating the kill probability based on an expansion of (l-z)^ 
in Jacobi polynomials rather than as a binomial series. The 
series is found to converge to the true value with less than 
n terms, which is very attractive for the calculation of 
salvo kill probability with a large n. 
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The author has tried to find alternative ways of calcu- 
lating the salvo kill probability which are effective even 
for the cases with large n's, but at present, the obtained 
results are not promising. The following propositions are 
given not as efficient algorithms for calculating the salvo 
kill probability, but as possible hints for developing more 
efficient ways of calculation. 

For the Elliptical Normal Salvo Model, we have a 
proposition similar to Proposition 2.5. We omit the proof, 
since it closely parallels that of Proposition 2.5. 

Proposition 7.2 

The salvo miss probability Q(n) of the Elliptical Normal 
Salvo Model characterized by (7.1), (7.2), and (7.3) with 

the weapon number n , 



Q(n) 



7T 



00 




00 



/ (1 




n 



-P 



X 



2 4.- 

s -p t 

y 



ds dt 



(7.8) 



satisfies the relationship 



Q(n) = I (5J) (l-y)^~^ q(k) , 

k=0 ^ 



(7.9X 



where 



q(k) = 



n/p 



x*^y 



7T 



/ / (1-e 



2 ^2 -p 

-s -t ^n^ Px 



2 .2 

s -P t 

^ ds dt 



(7.10) 



109 



In the case of Simple Salvo Model, q(k) had a simple 
recurrence relationship, but we have not succeeded in 
finding such a relationship in this case. 

Corresponding to the recurrence formula for the salvo 
kill probability presented in Chapter II, we have the next 
propositions; The series given in Proposition 7.3 is an 
infinite series, but for a small 6 /p it will converge rapidly. 

Proposition 7.3 

The salvo miss probability of the Elliptical Normal 
Salvo Model with n weapons is given by 



Q(n) = A -6^/p^ J • 2 J(n,2k) 

k=0 2^(k!) 



(7.11) 



where 



J(n,k) = e ^ t^ dt 



• 0 



with 



P 




(7.12) 



6 




110 



and 



y = 



a a 

— y 



^ ial + a^) (a^ + a^) 

X X Y Y 



Proof 

In (7.8) , let s = r cos 0, t = r sin 9, and (p +p )/2 = p, 

X y 

(p -p )/2 = 6. Then 
X y 



p s^ + p = r^(p + 6 cos 29) 

i^x ^y ^ 



and therefore Q(n) is 



Jp p 2 2 2 it r 2 

( M rr cos29^qj^ 



Q(n) = / (1 -ye"^ )"e"^^ {/ e‘ 



0 



= 2^PxPy (l-ye“^ )^e~^^ lQ(6r^) r dr , 



where Iq(z) is the modified Bessel f\mction of the first 
kind with order 0 , and has an expansion formula 



Io(z) 






k=0 



(7.13). 



Thus 
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oo 



Q(n) 



= /p p /p / (l-ye"^/P)^e"^I-(6t/p)dt 

X y 0 ^ 



^ k 2 
^ k=0 2 (k!) 



(2k) ! (5/p) 2k 



(2k) ! 



f°° -t/p> n -t.2k,. 

j (1 - ye ) e t dt . 

0 



□ 



Proposition 7.4 

The integral J (n,k) which has two integer arguments and 
two positive parameters y and p, 



J(n,k) 



1 f°° -t/p.n -t .k 

^ / (1 - ye ' "^) e t dt 

0 



(7.14) 



satisfies the following recurrence relation. 

n f “2t f 7 k ••• 

J(n,0) = ^ J(n-1,0) + ^ (l-y)"" , n = 1,2, 

J(0,k) = 1 , k = 1,2, . . . 



(7.15) 



(7.16). 



(7.17) 



112 



Proof 



00 

T/_ 1,N 1 ^ -t/pvIl-1 -t .k 

J(n,k) = j (1-ye e t 



) e t dt 



- y 



1 r i-[ ^“t/p n -1 -t/p-t 1 



kl 



dt 



The first term in the right-hand-side is J(n-l,k) . If 
the second term is integrated by parts, then 



J(n,k) 



T / T 1 \ 1 P /T ~t/p.n -t •} 

= J(n-l,k) - ^ (1-ye '^) e t 



1 P r t-i -t/p,n -t , .k-1 

FT n / ^ ® k t dt 



( 1 -ye"^/P)^e"^ t^ dt 

K. n Q 



Therefore, for k 7 ^ 0, we have 



J(n,k) = 



J(n-l,k) + ^ J(n,k-1) - ^ J(n,k) 



From this relation (7.15) is immediate. If k = 0, 



J(n,0) = J(n-1,0) + ^ (1 - y)^ - ^ J(n,0) , 



and therefore, (7.16) results. The relation C7.17) is obvious, 

□ 
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It is noted that all the terms in the expansion (7.11) 
are positive, and that the integral J(n,0) is the salvo miss 
probability of the Simple Salvo Model dealt with in Chapter 
II. Therefore we have 

Corollary 7.5 

A lower bound to the salvo miss probability of the 
Elliptic Normal Salvo Model characterized by (7.1), (7.2) 

and (7.3) with weapon number n is given by 



where Q(n) is the miss probability of the Simple Salvo Model 
by n weapons with 



Qg(n) = A - 6^P Q(n) 



(7.18) 



P 




6 




y 
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